IMPLEMENTATION AND TESTING OF A 
PROPOSED DIFFERENTIAL OMEGA SYSTEM 



Robert lawrence Vence 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 

IMPLEMENTATION AND TESTING OF A 
PROPOSED DIFFERENTIAL OMEGA SYSTEM 

by 

Robert Lawrence Vence, Jr. 

Thesis Advisor: C. E. Menneken 



June 1972 



Appnjoved {^ofi public. AeXai&e; diit^atcon usitunitzd. 



Implementation and Testing of a 
Proposed Differential Omega System 



by 



Robert Lawrence ^Vence, Jr. 

Lieutenant Commander, United States Coast Guard 
B.S., United States Coast Guard Academy, 1963 



Submitted in partial fulfillment of the 
requirements for the degree of 



ELECTRICAL ENGINEER 



from the 

NAVAL POSTGRADUATE SCHOOL 
June 1972 



ABSTRACT 



Omega- is a very low frequency (VLF) navigation system 
which will give world-wide coverage with eight stations when 
fully implemented. Using published skywave correction tables • 
accuracies of 1 to 2 n.mi. are attainable. Through the use 
of differential Omega, correction information can be dissem- 
inated to users in the vicinity of a monitor site. Differ- 
ential Omega accuracies are directly proportional to distance 
from the monitor site and are 0.26-0.5 n.mi. at 200 n.mi. 
from the monitor site. 

A system using U.S. Coast Guard radiobeacons as the dif- 
ferential information transmitter was proposed by Goodman 
and McKaughan. Their proposals were examined and improve- 
ments suggested. Two additional differential Omega systems 
being proposed by civilian contractors were also examined. 

The interface between an Omega receiver and the radio- 
beacon was built and tested. This device produced a CPE of 
25 yds. (0.16 cec) and a 95% circular error of 52 yds. (0.32 
cec) at the monitor site for LOPs crossing at 72® which is a 
marked improvement over accuracies that can be attained using 
skywave correction tables. ' 
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I. 



INTRODUCTION 



In mid-ocean, navigational accuracies of from 2-5 n.mi. 
are generally conceded to be adequate. As the coasts and 
harbor entrances are approached however, more precise navi- 
gation is desirable. This close-in region, generally called 
the coastal confluence region, extends to about 400 n.mi 
from the coast. Several systems are presently being used 
for navigation in the coastal confluence region and two of 
these are also useful for mid-ocean navigation. These 
systems are generally useful for ship positioning where high 
accuracy and reproducibility are required. Because mid- 
range frequencies are generally used, the range and accuracy 
of these systems will be limited during the evening hours by 
skywave effects. Lane ambiguity problems (to be described 
in detail later) are troublesome with phase comparison 
systems and require that the navigator know his initial 
checkpoint and keep an accurate track record. 

A. LORAN A/C 

Pulse time measurement systems are capable of fixing a 
ship's position by reception and identification of pulses 
received in a timed sequence from transmitting stations of 
known position. About 35% of the northern hemisphere is 
effectively covered by Loran A and 8% by Loran C [24]. 

Loran C equipment makes use of lower frequencies than Loran 
A and uses both pulse time and phase matching techniques. 
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Although the worldwide coverage is not as great for Loran C 

as it is for Loran A, the system is useful to greater ranges 

from the transmitter — 100-2300 n.mi. — and is capable of 

high accuracy. With good operating conditions Loran C 

accuracy can approach one foot per mile of range. Loran A 

/ 

has an estimated accuracy of 0.5 to 3.0 n.mi., depending on 
where the user lies in the hyperbolic net. Neither Loran A 
or Loran C suffers any ambiguity problems. 

B. RAYDIST "T” 

Raydist "T" is a hyperbolic phase comparison navigation 
system which requires four shore stations. Since this 
system uses phase comparison techniques for establishing 
ranges from the receiver to each shore station, lane 
ambiguities are a major problem. 

C. LORAC 

The Lorac (Long-Range Accuracy) system was developed for 
oil exploration activities. It uses a CW phase comparison 
method which requires three shore stations. The lane ambi- 
guity problem can be minimized by doxabling the number of 
transmitters. This system is used by both commercial and 
Government survey groups for near shore control work. Work- 
ing ranges are 100 - 150 n.mi. Ranging accuracies of 
1:50,000 are readily obtained with this system. 

Since Loran A was introduced to the navigator during 
WW II, the search for longer range, more accurate electronic 
navigation systems has been pursued. One of the proposed 
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Figure 1. Comparison of Navigation System:?. 
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systems, called Radux, used the arrival time of the signal 
in terms of the phase of the transmitted signal instead of 
the beginning of the pulse signal as Loran A does. As an 
outgrowth of Radux experiments and the shortcomings of pre- 
vious systems, a very-low frequency (VLF) navigation system 
called Omega was proposed. VLF was chosen because the Radux 
experiments showed that timing errors did not increase 
markedly with increasing distance and VLF showed extremely 
low attenuation rates [24]. The VLF band also shows very 
stable propagation characteristics that can be readily 
predicted. 

Despite the advantages of VLF transmission, there are 
several drawbacks, the most serious of which is the diurnal 
(daily) phase shift of the VLF signal due to changes in the 
ionosphere. Other effects caused primarily by the sun's 
emissions produce variations (anomalies) to the predicted 
diurnal phase values . 

These anomalous variations in the phase of the VLF signal 
are why the concept of differential Omega is of value. If 
some fixed site, referred to as the monitor site, receives 
the Omega signals and compares them to what they should be 
at the monitor site's geographic position, then the difference 
between the expected and the actual values can be disseminated 
by some means to Omega users in the immediate vicinity who 
may be subject to the same error. Upon receipt of this infor- 
mation the navigator may apply the corrections to his received 
Omega LOPs. The differential or relative Omega concept is 



I 



11 



not new but has been referred to before in "A Technical 
Evaluation of the Omega System" dated 18 December 1962 by 
NELC, San Diego and in many other reports since then. The 
use of differential Omega and Omega would allow one system 
to serve both the mid-ocean and coastal confluence navigator. 

Goodman's thesis [7] proposed the already established 
U.S. Coast Guard radiobeacon net as an inexpensive, readily 
available means of transmitting differential information 
along the United States' coastlines. Goodman's and 
McKaughan's [13] coding message formats have been examined 
in this thesis and improvements suggested. The design of a 
commutation circuit to implement the coding format for the 
differential Omega message was completed. A single line of 
position differential Omega correction device was built and 
tested first by itself and then using the CG radiobeacon as 
a transmission vehicle in a closed loop transmission system. 
The reasons for the anomalous phase variations were also 
examined. In addition other differential Omega systems 
currently being considered for implementation in the United 
States and France were investigated. 
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II. OMEGA 



A. SYSTEM DESCRIPTION 

Omega like Loran measures the difference in time of 
arrival of two synchronized transmissions and thereby 
arrives at a hyperbolic line of position representing an 
equal time difference. As mentioned previously Omega 
measures the time difference by comparing the phase of the 
arriving signal. Since the baseline — a great circle con- 
necting two Omega stations — is approximately 5000 - 6000 n.mi. 
long and the wavelength of 10.2 kHz radio waves is approxi- 
mately 16 n.mi., it is apparent that there are many wave- 
lengths of Omega signals and therefore many points of equal 
phase between the two stations. Line of position (LOP) 
ambiguities result because phase can only be electronically 
resolved into 360° segments within the receiver. The seg- 
ments are reduced from 16 to 8 n.mi. when navigating using 
the hyperbolic mode because, as a ship travels one-half a 
wavelength along the baseline between two stations, the 
phase of the received signal of the station toward which the 
ship is moving, compared to the receiver internal oscillator, 
will decrease 180°. At the same time the received phase 
comparison for the other station will increase 180° since 
the ship has moved one-half wavelength away from that station. 
The resultant phase-difference between the two stations is 
360° — this corresponds to an effective lane width for 10.2 kHz 
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on the baseline of 8 n.mi. Without additional information, 
the navigator cannot tell in which one of hundreds of 8 n.mi. 
lanes he lies. This problem is effectively eliminated by 
the transmission, of multiple frequencies as will be explained 
later. 

The extremely long baseline of VLF systems are especially 
effective on a spherical earth. For example, with a baseline 
subtending 60® of arc (approximately 4000 n.mi) the diver- 
gence between the hyperbolic LOPs is limited to a factor of 
two as one moves from the baseline instead of increasing 
infinitly as it does with a baseline negligibly long in 
comparison with the curvature of the Earth [24]. In addition 
to the decrease in divergence there is another important 
geometric factor which favors long baseline systems. Long 
baseline systems with optimum station placement insure that 
the crossing angles between LOPs will be large. 

The use of VLF transmissions allows more than the neces- 
sary three stations needed for hyperbolic navigation to be 

observed at any point. Because each of the eight Omega 

. . 3 

stations will serve about — of the surface of the Earth, 
there will usually be five or six stations providing useful 
signals [24] . Since any station can be measured with respect 
to any other, there will be either 10 or 15 LOPs available. 
The navigator is free to choose the three LOPs that give him 
the best signal or fix geometry. 

Unlike the Loran A system which requires the receipt of 
the master pulse at the slave station before the slave trans- 
mits, Omega transmissions are synchronized to a common 
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precise time standard, universal time. This change was made 

possible by the introduction in 1953 of atomic time standards 

accurate to one part in 10^°. Each Omega station has a 

bank of four cesium atomic standards which when intercompared 

. 1 2 

result in frequency accuracies better than one part in 10 
[21] . The standards are used both for carrier frequency 
synthesis and timing transmissions. 

The various Omega stations transmit in accordance with 
the format shown in Figure 2. The length of the transmis- 
sions can be seen to vary between 0.9 and 1.2 seconds. This 
transmission sequence is designed to enable the user to 
identify each individual station since the pattern surround- 
ing each station's transmission is unique. For example, to 
identify station D, the navigator would locate the 1.2 
second pulse preceeded and followed by a 1.1 sec pulse. 
Station D can thus be distinguished from station G which 
also transmits a 1.2 sec pulse. The beginning of station 
A's transmission is synchronized so that A transmits at 

0000 GMT each day and every ten seconds after that. There- 

» 

fore a user with a WWV or VT.WL receiver could manually 
synchronize his receiver commutator to match the transmis- 
sion sequence. Modern Omega receivers perform this 
synchronization automatically by measuring the length of 
the transmission pulses and comparing them to the known 
transmission pattern until station identification has been 
completed. 

Omega, as stated previously, differs from other hyper- 
bolic systems in the use of very-low frequencies (10-14 kHz) 
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Figure 2. Omega Signal Format. 



characterized by relatively stable propagation to very long 
distances (5000 to 8000 n.mi.). The 10.2 kHz signal provides 
the ship navigator with general navigation of the order of 
1 mile accuracy using the published skywave correction tables, 
but with 8 n.mi. ambiguities. This is not a limitation when 
tracking the 10.2 kHz signal continuously. If there is an 
interruption to the reception owing to a station failure, or 
failure of the receiver, the navigator can re-establish the 
correct 8 mile lane by other means. Dead reckoning would 
suffice during the Omega outage if the accumulated error 
remained less than 4 miles. Weather permitting, celestial 
fixes to this same order of accuracy would also define the 
correct Omega 8 mile lane. Failure of one station in the 
proposed 8 station network would enable the navigator to 
select another station without loss of lane identity. 

In order to reduce the dependency upon other navigation 
information to avoid lane ambiguity, additional frequencies 
are used. Use of 10.2 and 13.6 kHz provides a virtual 
difference frequency of 3.4 kHz. This lower frequency 
establishes a coarse lane with ambiguities spaced 24 n.mi. 
Although the precision of measurement at 3.4 kHz is reduced, 
it will suffice to identify the correct 8 mile lane within 
the 24 mile lane. For those navigators who know their posi- 
tion only to within ±36 miles, a third frequency, 11.33 kHz 
has been proposed. This would provide an unambiguous lane 
72 n.mi. wide when used' in conjunction with the 10.2 kHz 
signal. 



■ At the present time there are four Omega stations trans- 
mitting at reduced power. These are located in Norway, 
Trinidad, Hawaii and New York. The New York station will be 
replaced by one at La Moure, North Dakota before the summer 
of 1972. The following figure shows the anticipated goal 
for placing the present stations on the air at full power — 
10 kW — and the construction dates of the remaining four 

\ 

stations. 



10.2 kHz Trans- 
mission Segment 


Station 


On Air Date 
with 10 kW 


A 


Norway 


Late 1973 


B 


Trinidad 


Late 1973 


C 


Hawaii 


Late 1972 


D 


North Dakota 


Summer 1972 


E 


Reunion Is. 


Late 1973 


F 


Argentina 


Early 1974 


G 


Australia 


Early 1974 


H 


Japan 


Late 1972 



Figure 3. Omega Implementation Schedule [2] . 



B. VLF PROPAGATION 

Omega is a propagation limited system [28] . VLF prop- 
agation can best be regarded as being composed of various 
modes within the spherical waveguide formed between the 
earth and ionosphere. If several modes are present, modal 
interference may occur resulting in phase irregularities. 



18 



However, the principal mode may be expected to have the 
lowest attenuation rate and hence dominate at sufficiently 
large distances [23]. Because of this. Omega is considered 
a long range system and it is not recommended that any 
station closer than 600 n.mi. be used. The severity of 
modal interference will depend not only on separation from 
the transmitter but also propagation details such as the 
path, geomagnetic latitude, bearing and diurnal period. 

Since waves are propagated within the guide formed 
between the earth and the ionosphere, variation of modal 
parameters such as phase velocity may be expected as a 
result of changes in the ionosphere or ground. The obvious 
navigationally undesirable variation is diurnal phase change 
An uncorrected Omega LOP may vary as much as one complete 
lane from day to night. Ground-conductivity, latitude. 
Earth's magnetic field, path bearing, season and solar 
activity have also been demonstrated to have an effect on 
VLF phase. 

1. Skywave Correction Tables 

Diurnal and undesirable spatial variations are incor 
porated in skywave correction tables published by the Naval 
Oceanographic Office. All skywave correction tables pub- 
lished previous to 1970 were based on a global theory of 
Omega propagation incorporating theoretical and empirical 
physical principles where the relative contributions of the 
various effects are determined by regression analysis on 
millions of hours of data. The physical model has undergone 
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continual refinement at the Naval Electronics Laboratory 
Center for ten years [28], Since the fall of 1970 skywave 
corrections , limited areas have been receiving the benefit 
of a force fit wherein local prediction errors are determined 
by monitoring and then removed over whatever spatial extent 
may be justified by the statistics [10] . 

2 . Regular Diurnal Phase Change 

When dealing with changes of phase it is important 
to understand what is meant by phase decrease or phase 
increase. In all cases in this thesis, the term phase 
decrease will mean that the number of centicycles (cec) 
decreases.^ The term phase increase means the number of 
centicycles has increased. If one associates phase changes 
with changes of phase velocity brought about by changes in 
effective ionospheric reflecting height, a phase decrease 
corresponds to an increase in the phase velocity as a result 
of a decrease in height. 

The apparent height of the reflection layer during 
the day varies with the zenith angle of the sun . The for- 
mula is 

h = h(x=0) + h log secx (1) 

X s 

provided x < 85°. Here h(X=0) is the height for an over- 
head sun and h^ is the change in height of the ionosphere 

centicycle is a term which is used frequently in con- 
nection with Omega navigation signals . It is defined as one- 
hundredth of a full cycle of phase change at the frequency 
under consideration. Therefore one cec at 10.2 kHz equals a 
LOP displacement of approximately 480 feet on the hyperbolic 
system baseline. 
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at the region of reflection, has a mean value of 5.5 km 

but rises to about 7 km in October and April and falls to 
about 4 km in June and January. The value of h(X=0) is 
69 km [27] . 



long paths shows that the diurnal phase variation is trape- 
zoidal with almost constant values during periods of night- 
time or daytime propagation. A more or less steep phase 
increase and phase lag occurs at sunset and sunrise 
respectively. 

Diurnal phase change of VLF signals is commonly 
expressed in microseconds as 



where D is the path distance and and represent the 
night and day phase velocities of the first-order mode. 
This linear relationship is true only if no higher-order 
modes exist in the day or night earth-ionosphere waveguide 
at the receiver location [33]. Owing to the strong atten- 
uation of the higher order modes, the model of the single- 
mode propagation applies approximately to the propagation 
over long distances. 



Another expression used shows that the diurnal phase 



shift is almost directly proportional to the length of the 
path which is in daylight. This expression is 



The phase behavior of VLF signals propagated over 




( 2 ) 




( 3 ) 



21 













jf 



t 



I 




Figure 4. Typical Diurnal Phase Change (Haiku to 
Forestport, September 1968) [18]. 




.Figure 5.' Idealized Sunrise Transition on a Long 
East-West Path. Sunrise at the Trans- . 
mitter Occurs at Point A. The interval 
AB is approximately one-half hour. 

Phase shifts BC and DE are due to sudden 
changes in excitation factor. Point F 
shows sunrise at the receiver [18] . 
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where h = 80 km, is taken to be the mean height between day 
and night, Ah = change in height from day to night. 

At = change in time of arrival of signal and t = s/c where 
s = length of path in daylight under consideration and 
c = velocity of light [9]. 

This expression applies to the principal mode of 
propagation in a perfectly-conducting wave guide, and solu- 
tions using this expression are valid for distances greater 
than about 3000 km where the higher order modes may be 
neglected. 

a. Variation of Diurnal Phase Shift with Distance 
A final expression that expresses the change of 
phase of a received signal due to distance and the height 
changes in the ionosphere is 



A4> 




h 



(4) 



where X is the free space wavelength, h is the height of 
reflection, a is the radius of the earth. Ah the height of 
reflection changes, and d is the distance of the observer 
from the transmitter [5] . 

3. Other Anomalous Variations 
a. Solar Flares 

(1) Short Distance Effects . At short distances, 
up to 1000 km, the effects of solar flares on the phase of 
VLF signals are clearly visible. A flare effect is observed 
only during daylight hours and is associated with sudden 
ionospheric disturbances. These disturbances have a rapid 
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onset and slower recovery and last for times that vary from 
a few minutes up to a few hours. A phase decrease is 
observed at the same time as an HF radio fadeout and the 
optical sighting of the solar flare. It is found that the 
amplitude of the VLF signal does not change appreciably 
during the disturbance but the phase of the skywave decreases; 
hence the effect has been termed the "sudden phase anomaly" 
(SPA) . 

(2) Effects at Greater Distances . At oblique 

incidence, SPA's have been reported during the optical sight- 
ing of solar flares. The phase advance has an onset time of 

about 1 to 3 minutes , and recovery requires 30 minutes to 

2 or 3 hours [4]. The magnitude of the phase shift is 
apparently related to the increase in solar radiation, its 

angle of incidence at the ionosphere and to the length of 

the path over which a lowering of the apparent reflection 
height occurs. 

Sometimes solar flares occur which produce 
cosmic ray particles (solar protons) . The sudden phase 
anomalies observed during these events are unusually large. 
Such events are quite rare however. The rapid ionization of 
the lower region of the ionosphere as a result of and coinci- 
dent with the arrival of the cosmic-ray particles causes the 
phase of a VLF transmission to decrease at a rate apparently 
proportional to the increase in intensity of the solar cosmic 
radiation. 

(3) Origin and Energy Spectrum of a Solar Flare . 
During quiet conditions , the radiation of the quiet corona 
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and the transition region to the chromosphere exist perman- 
ently, whereas the emission of the hot regions occurs only 
during disturbed conditions. There is considerable evidence 
that when certain solar flares occur, the flux of energy 
arriving at the earth's atmosphere contains x-rays. These 

x-rays were produced in the coronal region at temperatures 

8 ... 

on the order of 10 degrees Kelvin. This is two orders of 
magnitude hotter than quiet sun conditions [4]. 

(4) Frequency of Flares Versus Disturbance to 
VLF Transmission . Just because a solar flare occurs, there 
is no direct probability that it will cause an SPA. For 
example from Ref. 1 it was shown that between late 1961 and 
mid 1964, 1846 optically detected solar flares were observed. 
Of these, 66 produced phase anomalies. The probability of 
occurrence increased approximately with flare area, but was 
less than probabilities for sudden ionospheric disturbances 
(SID) determined by other means. This indicates that VLF 
phase is not as sensitive as other SID indicators. 

The number of flares accompanied by SPA's 
varies from about 7 percent of all the observable flares in 
one month to none in another. For example, 104 flares were 
reported in June 1962; of these flares, 22 were of class 1 
or greater, but no VLF phase anomalies were observed [1] . 

(5) Teres tial Effects of a Solar Flare . See 

Figure 6 . 

b. Magnetic Activity 

Magnetic storms caused by solar flares also 
affect the phase of the downcoming skywave. The major 
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Figure 6. The Terrestrial Effects of a Solar Flare [15] 



difference is that these anomalies may occur by day or night 
and may last several days. These disturbances are associated 
with magnetic disturbances but not in a very definite way. 

The long distance VLF effects associated with magnetic 
activity appear to accompany isolated magnetic bays and 
great magnetic storms. During extreme magnetic storms the 
phase may change as much as 150° (42 cec) per minute. 

In addition to the above mentioned anomalies 
associated with great magnetic storms , smaller phase anom- 
alies are observed, during both day and night, which appeared 
to correlate with magnetic storms of less intensity. Some 
of these effects are found to occur about 3 days after the 
onset of the magnetic storm [27] . The magnetic variations 
generally occurred approximately 30 minutes earlier than the 
phase variations [4]. 

c. Meteor Showers 

In addition to solar and magnetic activity caus- 
ing phase anomalies, there is some evidence of an effect due 
to meteor showers. The anomalous effect in the diurnal phase 
pattern associated with meteor showers is shown in Ref. 4. 

The effect was observed during the hours of sunrise along a 
path from Great Britain to Boulder, Colorado when the phase 
is steadily decreasing (approximately 240° in 8 hours) . At 
the time of passage of the meteor shower, the steady phase 
decrease associated with sunrise is delayed until the meteor 
shower has passed, or until the entire path is in daylight, 
at which time the phase returns to the normal daytime value. 
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d. Atomic Device Explosion 

In order to understand what occurs to the phase 
of VLF signals during a high altitude nuclear explosion it 
is necessary to know how the D region is modified by high 
altitude nuclear bursts. There are several nuclear-burst 
processes which can greatly modify the D region of the 
ionosphere, including prompt-gamma-ray and thermal x-ray 
ionization, prompt-neutron decay, fission-beta decay and 
beta decay of the debris material [1] . 

(1) Prompt VLF Effects . The prompt VLF effects 
observed during nuclear testing occurred within seconds after 
the detonations. Since some of the paths considered are 
shielded by the earth from direct ionization radiation such 

as x-rays and yj^^ys created in the detonation, it is asserted 
that prompt neutrons from the burst decay into protons and 
beta particles, with a half-life of about 13 min; the betas 
then spiral down the geomagnetic field lines into the D 
region of the ionosphere where they enhance the ionization 
and cause increased absorption and decreased propagation 
times of VLF radio waves propagating in the earth-ionosphere 
wave guide. The geometry of the detonation point and VLF 
propagation path must be such that the geomagnetic field 
lines which intersect the propagation path are within line 
of sight of the burst [15] . 

(2) Delayed VLF Effects (1.5 to 5 minutes) . 

The VLF anomalies on paths where the geomagnetic lines are 
below the shadow cone which is tangent to the earth's 
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surface and whose apex is at the detonation point can be 
explained in terms of the eastward drift and energy spectrum 
of the electrons comprising the bomb-induced radiation belts 
which form in the geomagnetic meridian through the burst 
point. Those that are trapped will precess eastward at a 
rate proportional to their energy. A 1 MeV electron drifts 
at a rate of about 6-10 degrees of longitude per minute [1] . 

( 3 ) Delayed VLF Effects (15 to 40 minutes) '. 

VLF anomalies delayed by tens of minutes are believed to be 
caused by the eastward drift and energy spectrum of trapped 
radiation composed largely of neutron-decay electrons. It 
is possible that some of the field intensity anomalies 
observed may be attributed to trapped neutron decay electrons 
which drifted several times around the world [15] . 

(4) Phase Shifts of High Altitude Nuclear 
Explosions . High altitude nuclear explosions produce clearly 
identifiable phase shifts if the VLF path is located suffi- 
ciently close to the burst. The prompt neutron precipitation 
forms the burst signature — a sudden fast phase decrease, a 
slower decrease which lasts a minute or two, then a return 

to normal lasting tens of minutes, quite like that of a 
solar flare SPA. Figure 7 shows the magnitude of the phase 
shifts . 

e. Ionospheric Variations 

It has been shown earlier that changes in the 
height of the ionosphere alter the phase velocity of the 
signal in the earth-ionosphere waveguide. Therefore anything 
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Figure 7. Phase Shift of Omega Signals, Due to Three High Altitude Nuclear Bursts [16]. 



that affects the height of the D region, which acts as the 
upper surface of the earth-ionosphere waveguide, will* greatly 
affect the phase of the VLF signal. Since the occurrence of 
SID's follows the same 11 year occurrence cycle as do flares 
and sunspots, then their effect on propagation would be 
expected to follow the solar cycle. The season of the year, 
latitude of the observer and the position of the moon also 
have a direct effect on the phase variation of VLF signals 
[32]. 

There is an appreciable change in nighttime phase 
velocity with latitude. This change appears as an increase 
in the nighttime phase velocity with latitude for paths 
within the temperate zone. An apparent semidiurnal height 
variation of about 0.1 km appears to be related to the moon. 
The maximum height which occurs about two hours before lunar 
transit may be caused by a maximum rate of change in particle 
height which would be expected three hours before tidal 
maximum [32] . 

f. Polar Cap Absorption (PCA) 

Polar Cap Absorption is believed to be due to 
high energy solar protons in the 10 MeV range which penetrate 
deeply into the atmosphere over the polar cap and greatly 
enhance the ionization in the 50 to 80 km range [15]. The 
start of a PCA is from a fraction of an hour to several 
hours after a solar flare is observed. They usually last 
from 1 to 10 days with most lasting three days. Since the 
region of excitation is somewhat lower for these particles , 
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the general effect is not only a lowering of the ionosphere 
but also an increase in the .gradient of the reflection layer 
which usually results in an increase in attenuation rate 
[32] . Unusually large phase changes at VLF are observed 
during PCA so therefore it is fortunate that their occur- 
rence is rare. From October 1961 to October 1963 only seven 
PCA were observed on the polar path from Great Britain to 
College/ Alaska [1]. Since these are polar path events, 
they are of little concern to those who are using temperate 
latitude propagation paths. 

g. Ground Conductivity 

As the conductivity of the Earth's surface 
decreases, the phase of the reflection coefficient decreases, 
i.e., becomes a larger negative value. Watt has shown the 
manner in which the phase velocity changes with ground con- 
ductivities for day and night paths. His conclusions were 
that the change in velocity is appreciably greater for the 
lower frequencies and also that the nighttime changes in 
relative phase velocity are appreciably less than those that 
occur during the day [32]. 

Swanson compared the variation of phase velocity 
with electrical ground conductivity for the daytime ionosphere 
and also the velocity variation of 10.2 and 13.6 kHz signals 
with ground conductivity [17] . He concluded that the varia- 
tion of velocity at 10.2 kHz was relatively minor until 
conductivities become quite low at which point variations 
become more significant. In the case of 10.2 kHz the phase 
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velocity decreased monotonically as the conductivity 
decreased. In the case of 13.6 kHz, however it appeared 
that a rounding or change may be occurring in the region 
0.006 to 0.1 millimhos per meter. The large anomaly 
centered near 20 kHz apparently shifts toward lower fre- 
quencies as conductivity is decreased. The effects are 
therefore significant and cannot be anticipated with any 
certainty without direct calculations for the conductivity 
range involved [17] . Since skywave corrections are now 
being force fit to meet local conditions, the effect of 
ground conductivity changes will become minimal, 
h. Geomagnetic Path Orientation 

The Earth's magnetic field has an appreciable 
effect upon the motion of the electrons in the lower iono- 
sphere and this change in motion can affect the relative 
reflection characteristics for VLF waves. There is not 
only a significant change in the amplitude of the reflection 
coefficient with direction of propagation, but there is also 
an appreciable change in the phase of the ionospheric 
reflection coefficient as the direction of propagation changes 
from east to west. The general tendency is for the phase to 
decrease as the direction of propagation changes from west- 
east to east-west, i.e., to the west propagation has less 
phase lag than to the east [32]. This can be explained by 
the radio wave interacting with the Earth's magnetic field 
resulting in an aiding effect east to west and hindering 
effect west to east. The attenuation rates are also 
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different for west-east and east-west transmissions at VLF 
with to the west being less. Since the characteristics of 
the ionospheric reflection coefficients are a function of 
the sharpness of the ionospheric boundary as well as fre- 
quency, the directional properties are different for day and 
night conditions. Nonreciprocal effects also change through- 
out the frequency range of the VLF band. Since the charac- 
teristics of the ionosphere can change throughout the solar 
cycle, directional effects also change throughout the solar 
cycle. The two primary influences which are important in 
this regard are the change in collisional frequency with 
solar activity and the change in galactic cosmic ray 
intensity [32] . 

The NELC computer program used for obtaining the 
skywave corrections and curves which approximate the complete 
wave solution, uses an expression of the form: 

b^sin 9 + B cos 20 + K1 = bsin 0 - 2Bsin^0 + K2 (5) 

where 0 = bearing angle with respect to magnetic field, 
bj^, B, K1 and K2 = constants. 

Swanson also compares the experimental curve 
with that obtained from the NELC program. There is quite 
good agreement between them. However, the NELC computer 
results of the complete wave solution are applicable only 
for conditions similar to those between Hawaii and San Diego. 
The experimental results are based on the synthesis of data 
involving not only a wide range of bearing angle, but also 
various values of the total field and various dip angles [17] . 
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i. Propagation Over Land-Sea Boundary 

When the propagation path is mixed, i.e., con- 
tains areas of appreciably different conductivity, determin- 
ing the phase and amplitude of a signal across the boundary 
becomes more difficult. When the distance separating the 
transmitter and the boundary discontinuity is relatively 
small, it can be shown that the transition from the form of 
variation with distance characteristic of the first medium 
to that of the second medium is quite rapid. On the other 
hand, when the distance from a transmitter is large, the 

I 

transition takes place over a much longer distance. At VLF 

the phase is very sensitive to the difference in conductivi- 

_2 

ties between land (a = 10 mho/m) and sea (a = 4 mho/m) . 

For example at 1000 km an all land path experiences an 
angular shift of -82° whereas an all sea path has only a 
-33° angular shift [32] . 

j. Fluctuations Due to Whistler Mode Signals. 
Whistler mode signals are signals radiated into 

the ionosphere near the transmitter and guided along the 
Earth's magnetic field lines back to the Earth in the 
vicinity of the magnetic conjugate of the transmitter. The 
magnetic conjugate is the point where the magnetic field 
line reenters the Earth's crust. A receiver near the conju- 
gate point may thus receive the normally propagated signal 
together with the whistler mode signal. Since the whistler 
mode signal varies with time, the resultant signal will 
exhibit, in general, both phase and amplitude variations. 
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Some observations have been made in New Zealand 
of the phase stability of whistler mode signals at VLF from 
NPG/ Seattle. The magnetic conjugate of NPG lies about 
3000 km SE of Lower Hutt, New Zealand and its direction is 
nearly perpendicular to the great circle path from NPG. 

Thus a pair of perpendicular loops were used to receive the 
direct and whistler signals separately. The observations 
indicated that the whistler signals tend to be shifted in 
frequency by several cycles per minute, the frequency shift 
being greatest (10 Hz/min) near midnight. The signals were 
observed when the source and conjugate were dark, but not 
when both were in daylight. The amplitude of the whistler 
mode signals was about 10 percent of the direct signals [4] . 

It is thus apparent that phase observations on 
NPG at Lower Hutt using an omnidirectional antenna would be 
characterized by phase fluctuations of the order of ±6® 
occurring several times per minute. Such variations might 
be important in a phase recording system having a fast 
response, but would be removed with phase recording systems 
with longer time constants such as the Tracer 599R with a 
time constant of 70 seconds. 

4 . Model for Phase of Received Signal 

As delineated earlier in this paper, the phase velo- 
city at VLF and hence 10.2 kHz is dependent upon numerous 
parameters. The parameters are again: 

a. Propagation over sea (distance) 

b. Propagation over land (distance) 
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c. Ionospheric height caused by 

(1) Day or night 

(2) Sun's zenuth angle 

(3) Solar activity 

d. Magnetic abnormalities 

e. Earth's magnetic field 

f. Earth conductivity 

g. Latitude effects 

h. Excitation factors 

i. Season of the year 

E. R. Swanson of NELC has developed a prediction model 
and a computer program to determine the parameter coeffi- 
cients for that model [15]. 

In general, four regions in the electromagnetic 
field may be considered at 10.2 kHz: near field — 0 to 30 km; 

medium field — 30 to 1000 km; far field — 1000 km to the 
boundary of the antipode region; and the antipodal region 
approximately 1000 km wide. The antipode is the point on 
the surface of the Earth diametrically opposite the trans- 
mitter. The near and medium fields are contaminated by 
high-order modes, skywaves, ground-wave interference, and 
the induction field of the transmitting antenna. The anti- 
podal region will be excluded because of interference between 
the same signal coming from several directions simultaneously. 
This leaves only the fer field which can be considered on a 
linear basis, permitting simple summation of propagation 
effects. Thus, the propagation can be a function of several 
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parameters such as time, distance, and others mentioned 

above. Therefore the equation for a plane wave would be 

2TTi(ft-k (d,t)D) 

E = |e^ (d) I e (6) 

a 

and the phase of a signal in cycles at some point is: 

4 . = ft-k^(d,t)D (7) 

However, k^ is generally independent or, at least, 
nearly independent, of time so 

k^(d,t) = k^(d) 

and 

^ 3k (d)D 
c 

3D 

Since, in the far field, the parameter coefficients 

were assumed to be linear the coefficient k (d) can be con- 

c 

sidered to be the sum of several terms and a correction 
factor: 

kc(d) = K^D + + 4>^,(d) (9) 

and may be determined by linear regression and 
is a constant and lost upon differentiation of equation (9). 
Furthermore, a phase designation also independent of time 
can be substituted for k^(d) in equation (9), resulting in 

P = K^D + K^, cycles (10) 

is an excitation factor which merely indicates the phase 
lead or lag of the field relative to the antenna. Equation 
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(10) permits the calculation of phase at any point but 
assumes like conditions for the entire path. However, it 
is very unlikely that a uniform path will ever be encountered 
in actuality. So, equation (10) must be expanded. Since 
the far field only is considered in Omega navigation, linear 
propagation paths are assumed and of equation (10) 
becomes : 

n 

K, = Z K. (11) 

^ i=l ^ 

and for a land-sea path: 

P = K D + K,D, + K (12) 

s s 1 1 o 

where 

K = Parameter coefficient over sea water 
s 

= Parameter coefficient over land 
D = Normalized distance over sea, referenced to 10.2 kHz 
= Normalized distance over land, referenced to 10.2 kHz 
K = Excitation factor 

o , 

Equation (12) requires expansion to be really useful, 
since it includes only expression for land and sea propaga- 
tion and excitation. It neglects important parameters such 
as ground conductivity, sun's zenith angle. Earth's magnetic 
field and others. However it can be expanded to become gen- 
erally applicable; 

P = K + + ...• + K B (13) 

0 112 2 n n 

where defines respectively the parameter coefficients 
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and the parameters. When the individual parameter's expres- 
sions and their coefficients are determined, equation (13) 
can accomodate them and will result in increased accuracy of 
position. 

NELC, San Diego has determined a number of the 
coefficients experimentally and is using these for tabula- 
tion of skywave corrections. 

The form of equation (13) presently used is 

P = K,D, +KD +K + + K'B' +K^B^ (14) 

11 ss o 55 55 66 

where K^,D^,K^, and were previously defined and: 

= Coefficient of magnetic orientation 

Bg = Parameter of magnetic orientation /cos<j)sin6dD 
= Additional coefficient of magnetic orientation 

Bi = Additional parameter of magnetic orientation, 
/cos(})sin^ 6 dD 

Kg = Coefficient of Latitude 

Bg = Parameter of latitude, /cos(|)COS (dip) dD 

where 6 = magnetic bearing and (J> = geomagnetic latitude in 
degrees . 

The values of the coefficients can be found in table 
3.1 of Ref. 15. 

Equation (14) has been derived in such a manner that 
it can readily include any additional parameters defined at 
some future date. Accuracy of the phase model, and conse- 
quently the fi-X position, will improve as pertinent parameters 
are defined and included. Omega system operation will 
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furnish the data from which these additional parameters can 
be determined. 
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III. DIFFERENTIAL OMEGA 



A. GENERAL CONCEPT 

The effectiveness of Omega is predicated on the predict- 
ability of a VLF phase measurement made at any time and 
place. Environmental influences active over given areas 
for various periods of time contribute to anomalous phase 
variations. Use of published skywave corrections [31], which 
are based on a time and location dependent propagation model, 
eliminate a large portion of the anomalous variations. How- 
ever, certain residual phase errors remain and result in 
Omega position inaccuracies on the order of one n.mi. 

Because of the space and time characteristics of both 
Omega propagation and the skywave correction model, it is 
natural to expect that any remaining phase errors will to 
some extent be correlated in distance and time rather 
than simply being a random process. It is seen from Figure 8 
that a random, unpredicted ionospheric disturbance in area A 
will cause similar or correlated phase variations at points 
PI and P2 because the paths are effectively parallel and 
closely spaced. Since the distances involved with Omega are 
quite large (propagation paths in excess of five thousand 
miles are possible) , the separation between PI and P2 could 
be appreciable while remaining errors would still be 
correlated. 
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Figure 8. 



Effect of Ionospheric Disturbance. 
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If PI is a fixed geographic site whose Omega coordinates 
are accurately determined, then this site — called the monitor 
site — can compare the received Omega signal to the correct 
value associated with an undisturbed path and obtain a 
"difference" between actual and correct values. This is the 
meaning of the term differential Omega. If the monitor site 
has the capability of transmitting this differential infor- 
mation to a user in the differential operating area, then a 
complete differential Omega system is formed (see Figure 9). 
The size of the differential area has not yet been experi- 
mentally determined — however it can be said that the area 
would be no larger than the area in which the use of the 
differential correction as transmitted by the monitor site 
provided better results than the use of Omega and skywave 
correction tables. Once the skywave correction results were 
equal or better than differential Omega then only skywave 
corrections would be used. Different investigators have 
proposed various maximum separation distances for PI and P2 
but it is generally agreed that 250 n.mi. is a good approxi- 
mate value [19,20,25,30]. 

B. DETERMINATION OF SPATIAL CORRELATION OF OMEGA SIGNALS 

From the intuitive notion of differential Omega investi- 
gators were led to an examination of Omega phase error data 
to determine an- approximate model for the spatial correla- 
tion. Separations between locations where extended data was 
available were too large to adequately analyze differential 
Omega in a straightforward line-for-line approach. Instead 



44 




45 



Figure 9. A Differential Omega System. 
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Kasper [6,11,12] used statistical analysis techniques and 
autocorrelation functions to describe Omega error spatial 
correlation. The function used was; 

, -|4x|/1500 

R„(Ax) = e (15) 

where Rj^(Ax) = E [6(J> (x) 6([) (x+Ax) ] and E(*) denotes expected 
value. The users' location is defined as x and the monitor 
site as x^. The uncorrected error at x is 6(|)(x); the cor- 
rection broadcast by the monitor site is the phase error at 
x^, 6(|)(x^). Thus the residual error after correction is 

A<})(x) = 6({)(x) - 6(}){x^) (16) 



RMS error is found by squaring 

A(j)^(x) = 6({)^{x) - 26<}) (x) 6(}) (x^) + 6ct)^(x^) (17) 

taking expected values. 



E[A(J)^(x)] = E[6(})^(x)] - 2E[6(}) (x)6(|) (x^) ] 



+ E[6<\> (x^)] 



5 . -|AxI/1500 , 

= 4 - ^4 ^ ^ 4 



and square roots 



2-2e 



-| Ax|/1500\'^ 



(18) 

(19) 

( 20 ) 



where Ax is the separation distance between the monitor and 
the user in n.mi. RMS differential Omega error is shown as 
a function of Ax in Figure 10 (upper curve) . A value for 
Of. of 1 n.mi. is used. Within 200 n.mi. of the monitor 
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Equation (20) 
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Figure 10. RMS Differential Omega Error. 



site, RMS errors are less than 0.5 n.mi. Receiver errors 
are not included in the above analysis. 

Another statistical model for Omega error has been con- 
sidered [12] ; its autocorrelation function is 



R(Ax) 




I Ax|/698 
e 



1 + 



lAxI 

698 



( 21 ) 



RMS differential Omega error for this model is obtained the 
same way as before resulting in equation (22) 

-IaxI , 

1 + 4^1 ( 22 ) 






2-2e 



698 



This function is also shown in Figure 10 (lower curve) . 
Within 200 n.mi. of the monitor, RMS differential Omega 
errors are less than 0.26 n.mi. if a 1 n.mi. RMS unaided 
Omega error is assumed. 

It is likely that the actual error behavior lies between 
these two results. Gossard and Norgard [19] stated that 
the accuracy of differential Omega is inversely proportional 
to receiver separation and will be 0.2 statute miles at a 
separation of 100 miles. These investigators arrived at 
their conclusion in the following way. The bearing error, 
A6, can easily be related to the phase deviation, A(|), as 
shown in Figure 11. It can be seen that 



A(}) = -^ 6 = ^ dsin AS = ^ dA6 



(23) 



for small bearing errors. It follows that 





2tt 
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(24) 



Transmitter 




Figure 11. Relation of Bearing Angle Error 
to Phase Deviation. 
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For small fluctuations A(J> is proportional to d so the 



slope — — ^ = constant. Eventually, the separation of 

d 

receivers (1) and (2) will become so great that the phase 
fluctuations will be uncorrelated and 

= /T ‘ constant (25) 

By taking the slope of the RMS phase difference vs. receiver 
separation they determined the differential omega fix error. 

The above results are dependent only on the angle between 
the two receivers and the distance between them. It appears 
that this result will give the phase difference due to dis- 
tance and angle between the two locations but will not give 
any information about the effects on phase of SIDs, PCAs or 
magnetic storms along tbe propagation path. 

C. PROPOSED DIFFERENTIAL OMEGA SYSTEMS 
1 . French System by Sercel Company 

The Sercel Company has been experimenting with a 
differential Omega system which showed great promise in tests 
made in late 1971 over distances of from 5 to 300 km from the 
monitor station. The results will be discussed later. The 
Sercel system requires a fixed monitor site from which to 
obtain and transmit the variations in the signals received 
from each Omega transmitter. There are as many correction 
values to be transmitted to the mobile receivers as there 
are useable Omega signals in the area. The phase determined 
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by the receiver at the monitor site for each signal is with 
respect to the phase of its local oscillator. The local 
oscillator in the monitor and mobile receivers must be very 
stable so that its frequency doesn't change appreciably from 
the time that one signal is received until the second signal 
is received (as much as five seconds) . If it did change 
then the phase reference would be different for the two 
signals and some LOP error would result. Since at the mobile 
receiver all the received Omega signals are with 
respect to the same local oscillator/ when two stations' 
phase signals are subtracted from each other to generate a 
hyperbolic LOP, the common local oscillator reference signal 
is eliminated. Phase correction signals at the monitor 
site for the respective stations are continuous 1 kHz sine 
waves. The value of the phase correction to be transmitted 
is in the form of a phase difference with respect to the 
1 kHz internal reference, the latter in turn being phase 
coherent with the local oscillator. In effect, the analog 
transmission contains a total of five signals; four phase 
signals plus the reference signal. In order to minimize 
the transmission bandwidth all the signals are first reduced 
to a conveniently low (20 Hz) frequency, which constitutes 
a satisfactory compromise between the full transmission 
width and the effect of transmission defects likely to 
result from transmission speed and signal fading [14]. All 
five signals are then time multiplexed so as to correspond 
with the basic Omega format. This multiplexed information 
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Figure 12. Sercel Differential Omega Transmitting Station. 
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Figure 13. Sercel Differential Omega Receiver Station. 
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is then routed to a 100-watt transmitter in the 2 MHz range 
where it phase modulates the carrier. 

The shipboard installation includes a correction 
receiver which accepts and handles the phase correction 
signals. This receiver is connected via a single cable with 
the standard Sercel Omega receiver. The correction receiver 
contains a double superheterodyne section, the second 
frequency converter being driven by a small frequency 
synthesizer. This offers a choice of four frequencies 
above and four below the reception frequency, the interval 
between any two frequencies being 160 Hz. The second channel 
has a bandwidth of 50 Hz adapted to the frequency spectrum 
of one transmitter. The receiving system is thus capable 
of receiving the signals from eight correction transmitters 
which may be installed in a chain up and down the coast. 

Selection of a station may be achieved by simply 
selecting the frequency of the second local oscillator in 
the receiver. A discriminator and a filter placed after 
the second frequency converter is used to separate the 20 Hz 
signals. The timing signals obtained from the ship's Omega 
receiver are used by the correction receiver as demultiplex- 
ing commands for the 20 Hz signals. The four 20 Hz signals 
are demultiplexed and demodulated according to two quadrantal 
components by a 20 Hz signal generated locally in the cor- 
rection receiver. The demodulated signals are routed through 
narrow-band filters and the resulting voltages appear across 
four pairs of leads (one pair per station) . These voltages 
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are equal to the sines and cosines of the angles correspond- 
ing to the corrections. 

In conventional Omega operation received signals are 
applied to the receiver's meters which show the phase dif- 
ferences between the selected stations. In differential 
Omega navigation the four signals which carry the phase 
information of the Omega station being received are taken 
from the Omega receiver and routed over a cable to the 
correction receiver which contains four electronic phase 
shifters. Each of these is driven by voltages obtained from 
the correction receiver equal to the sines and cosines of 
the correction angle. Each of the signals taken from the 
Omega receiver goes through the corresponding phase shifter 
where it is subjected to a phase shift equal to the correc- 
tion angle prior to being returned to the Omega receiver 
meter input. By this method corrections are received and 
decoded by the correction receiver and automatically applied 
to the ship's Omega receiver. Figure 14 shows the results 
obtained from the tests held in late 1971 for two hyperbolic 
LOPS/ Norway-New York (A-D) and Trinidad-New York (B-D) . 
oA-D and oB-D are the standard deviation of the LOPs from 

the correct values, o, is the standard deviation of the fix 

a 

errors from the known fix position. 

2 . Micro-Omega 

Micro-Omega is the name of a differential Omega 
system developed by Teledyne Hastings-Raydist/ Hampton, Va. 
The operation of the Micro-Omega system is quite similar to 
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ad, meters with 
lane A-D = 23 km 
lane B-D = 44 km 
crossing angle 80° 
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Figure 14. Test Results of Sercel Diff. Omega System [14], 
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the Sercel system although it is not as adaptable as the 
Sercel system and uses a different modulation frequency and 
technique. In the Micro-Omega system. Omega signals 
received at the monitor station are converted to low fre- 
quency (300 Hz) audio tones, phase locked to one of the 
incoming signals. Each of the audio tones is phase shifted 
a fixed predetermined amount, which has the effect of sub- 
tracting out the known Omega readings. The audio tones are 
then retransmitted as single sidebands on an HF carrier. 
Omega signals received on the vessel are processed similarly 
to obtain a sequence of audio tones of the same frequency. 

A vessel using Micro-Omega receives the single side 
band transmissions. The latter are processed to recover 
the audio tones. Thus, a pair of simultaneous audio tones 
at a single frequency is obtained on the vessel for each 
Omega transmission. The audio tones resulting from a single 
transmission are applied to a phasemeter. Unique phase 
angles are thereby derived for each station. The phase 
angles with respect to the local oscillator for two stations 
are then fed to a subtractor circuit which gives the phase 
difference and generates the hyperbolic LOP on which the 
vessel lies. 

The analysis below from Reference 8 derives the 
phase relationships which exist within the Micro-Omega 
system. The formulas derived only apply to one line-of- 
position but are equally applicable for as many Omega sta- 
tions received at the monitor and user locations. 
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Figure 15. Micro-Omega Transmitting Station. 




Figure 16. Micro-Omega Shipboard Receiver 
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Figure 17. Micro-Omega System Geometry. 



The following definitions are used; 



0T = Instantaneous phase of signal A at 
transmitter A 

0Tj^ = Instantaneous phase of signal B at 
transmitter B 

0T^R^ = Instantaneous phase of signal A at vessel 

0Tj^R^ = Instantaneous phase of signal B at vessel 

0T R = Instantaneous phase of signal A at Micro- 

^ Omega Sta. 

0T, R = Instantaneous phase of signal B at Micro- 

^ Omega Sta. 

Dl,D2 ,D3 ,D4 ,D5 = Distances between stations shown in 

Fig. 17 

c = Propagation velocity of radio waves in air 

t = Time from some arbitrary initial time, 
t is eliminated later in the algebraic 
manipulations 

el,e2,e3,e4 = Distance error equivalent of skywave along 
paths 1,2/3, and 4 

N = Multiplying factor in the phase locked 
loop local oscillator 

We have 0T = 0T, = 27rft 

a b 

0T R = 2TTf t - 2TTf t - (Dl+el) (26) 

3. V C C C 

0TaRm = 2TTft - (D2+e2) (27) 

At the vessel the local oscillator is synchronized to the A 
Omega signal and assumes a frequency and phase relationship 
as follows due to the multiplier N. 

Local oscillator ^ N / . 2tf ( 23 ) 

phase on vessel N-1 I cl 
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At the Micro-Omega Station the same synchronization on the 
A Omega station is accomplished (see Fig. 15) . 



Local oscillator phase N 2-nf 

u- (27rft (D2+e2)) 

at Micro-Omega Station N-1 c 



( 29 ) 



and the following signal is fed to the transmitter 

^ (2arft - ^ (D2+e2)) (30) 

Both oscillators remain stable through B transmission when 
Omega B transmits. For the B transmission, 

0T, R = 27rft - (D3+e3) (31) 

D V C 

0T, = 2?rft - — (D4+e4) (32) 

Dm c 

When 0T, R is received by the ship and 0T, R is received by 
the monitor site, they are heterodyned with their local 
oscillators. The difference signal between the received 
signal and the local oscillator is fed through the commutator 
(see Fig. 16) to the phasemeter on the vessel. The corres- 
ponding signal at the Micro-Omega station modulates the 
Micro-Omega carrier. 

The IF output = Local oscillator - Incoming signal 
On the vessel direct from B 

IF output = (2»'ft - (Dl+el) ) - 2yft 

+ (D3+e3) (33) 
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IF output = 27rft- (Dl+el)) 

+ (D3+e3) (34) 

c 

On Micro-Omega Station direct from B 

IF output = (2irft - (D2+e2) ) 

- 2irft + (D4+e4) (35) 



IF output = (2Trft) - (D2+e2)) 



+ (D4+e4) (36) 

The latter signal, equation (36) , is transmitted to the ves- 
sel as modulation from the Micro-Omega Station. The vessel 
receives a slightly delayed signal at the lower frequency 

(jj^ (2iTft)) (37) 

When B is "ON" the vessel receives from the Micro-Omega 
station 



1 

N-1 



(2irft) 



1 .2irfD5 
N-1 ^ c 



N , 2iTf 
N-1 ^ c 



(D2+e2) ) 



+ ^ (D4+e4) (38) 

When A is "ON" the vessel receives from the Micro-Omega 
station 



^ (2iTft - ^ (D2+e2)) 



1 .2irfD5. 
N-1 ^ c ^ 



(39) 



61 



Comparing both A signals (equations (26) and (39)) in a 
phasemeter gives 



^ (^ (Dl+el-D2-e2-D5) ) 

N-1 C 

Comparing both B signals (equations (31) and (38) ) in 
phasemeter gives 

N ^2TTf (Dl+el-D2-e2) ) + (D4+e4-D3-e3) 

W J- c c 



( 40 ) 



a 



1 /2TTfD5v 

Subtracting A phasemeter from B we have 

(D4+e4-D3-e3) + (Dl+el-D2-e2 ) (42) 

c c 

Now since e4 = e3 and el = e2 and D4 and D2 are constants 

Phasemeter = (D1-D3) (43) 

c . 

This phasemeter output is equal to the correct hyperbolic 
line of position with no skywave error included. 

Published test results of the Micro-Omega system are 
quite limited. The data obtained during an eighteen month 
evaluation program yielded an RMS system error of less than 
one centi cycle (480 ft on Omega baseline) under favorable 
signal-to noise conditions. Reference 8 stated that these 
conditions prevailed more than ninety-five percent of the 
time for one line (Trinidad-New York) and only twenty per- 
cent of the time for the other available line (Hawaii-New 
York) . There was no data presented as to position errors 
versus distance from the monitor site. 
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IV. DIFFERENTIAL OMEGA SYSTEM USING COAST GUARD RADIOBEACON 



A. GENERAL DESCRIPTION 

Goodman [7] first proposed using the present system of 
U.S. Coast Guard radiobeacons, Type T-747A/FRN, as the 
transmission vehicle for the differential correction message 
along the coasts and harbors of the United States (see 
Fig. 18) . The correction information would be in the form 
of audio tones that are directly related to the Omega cor- 
rections for specific hyperbolic LOPs. A total of four LOP 
corrections and five calibration tones would be transmitted 
in a one minute interval once every five minutes. The 
message format will be described in detail later but in 
general a 600 Hz tone would be equivalent to no LOP correc- 
tion and 500 Hz excursions above and below this value would 
correspond to +50 cec LOP corrections. This audio tone would 
be used to AM modulate the radiobeacon carrier (285-325 kHz) . 

Instead of using the standard AM modulation technique 
where the carrier's amplitude is modulated with the signal 
of varying frequency, Goodman proposed using the dual carrier 
technique presently being used on many Coast Guard radio- 
beacons. This modulation method uses two RF carriers sep- 
arated in frequency by the audio frequency corresponding to 
the audio tone desired. When received these two carriers 
produce the same effect in the receiver as conventional AM 
but require only half the bandwidth — thus conserving frequency 
spectrum. 
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Figure 18. Goodman's Differential Omega 
Transmission System. 
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The use of the radiobeacon and the differential Omega 
message proposed by Goodman were predicated on the assumption 

I 

that a small vessel could use the system as well as the large 
vessel having computer oriented navigation systems . In order 
to use the proposed system, the shipboard user must have a 
low frequency (LF) receiver aboard as well as an Omega 
receiver (see Fig. 19). For the small users (commercial 
fishermen, coastal freighters, etc.) the audio output of the 
LF receiver is fed into a simple frequency meter that has a 
range of 100 to 1100 Hz. The frequent calibration tones used 
by Goodman were intended primarily for the navigator using 
an analog frequency meter. The face of the meter could be 
calibrated directly into corrections according to the con- 
version scale described below without having the navigator 
manually convert the frequency into corrections. These 
corrections can then be applied to the received Omega 
readings . 

As the sophistication of the shipboard navigation system 
increases, a digital meter can be used instead of the analog 
meter used above. The digital meter would be preferred 
because of its higher accuracy and ease of reading in a 
shipboard environment. A standard digital frequency meter 
is easily modified to read directly in centicycles instead 
of frequency. For Goodman's proposal all this would entail 
is subtracting 600 from the frequency received and dividing 
the result by ten. 
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CORRECTED OMEGA LOPS 
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For automatic operation/ the output of the LF 
receiver could be fed to a small special purpose computer 
which may include a frequency counter. The output of the 
Omega receiver can also be routed to the computer and the 
corrections applied to the Omega readings and displayed on 
the computer console or rerouted back to the Omega receiver 
and displayed on it. The special purpose computer would 
have to have an input register which identified the LOPs to 
which the correction message (described below) applied and 
what LOPs were being fed into it from the Omega receiver. 
The initial calibration tone could be used as an indication 
of the start of the message and the LOP identifiers as an 
indication to the computer as to which LOP it was to apply 
the correction. The correction could be determined by the 
same method described above in connection with’ the digital 
frequency meter. The computer operations can also be 
handled by a sub-routine in a large, multi-purpose computer 
already aboard instead of a special purpose computer. The 
important point that should be reemphasized however is the 
automatic system can use the same correction message as the 
manual user and thus the proposed system does not discrim- 
inate against the small, unsophisticated user. 

B. CHANGES TO GOODMAN'S PROPOSAL 
1. Transmitter 

The T-747-A/FRN . radiobeacon transmitter can be used 
in a variety of ways without additional circuitry: 

a. Keyed carrier, no modulation 

b. Keyed carrier and tone modulation 
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c. Continuous carrier with keyed tone modulation 

d. Dual carrier, both keyed simultaneously 

e. Dual carrier, either keyed separately 

Goodman also investigated the following ways of modulating 
the radiobeacon — all of which would require transmitter 
modifications : 

a. PAM — Pulse Amplitude Modulation 

b. PWM — Pulse Width Modulation 

c. PPM — Pulse Position Modulation 

d. FM — Frequency modulation 

e. Amplitude modulation with a modulation signal 
of varying frequency 

f. Carrier separation Modulation (similar to dual 
carrier modulation above except that one of the 
carrier's frequencies varies). 

Goodman decided that the carrier separation modula- 
tion technique was best since it only required changing the 
second RF carrier oscillator in the transmitter to a voltage 
controlled oscillator (VCO) and secondly used only half the 
frequency spectrum of conventional AM modulation. 

One of Goodman's major goals was the reduction of 
the number of modifications required to the radiobeacon 
transmitter. In this thesis amplitude modulation with the 
modulating signal of varying frequency is suggested since 
there are no wiring changes required to the transmitter under 
the arrangement described below. 

To implement the system proposed in this thesis all 
that is required is the removal of the temporary jumper con- 
nection between JlOl and J102 or J102 and J103 on the 
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Figure 20. Proposed Radiobeacon Differential Omega System. 
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transmitter osci,llator chassis. The connection of the Omega 
correction device (described below) to J102 is the only 
other step necessary for operation of the transmitter as a 
differential Omega transmitter. It should be noted that the 
radiobeacon can still be used by the non-Omega user for 
direction finding purposes since the radiobeacon normally 
is AM modulated with a 500 or 1020 Hz signal. The addition 
of differential Omega information merely makes the modula- 
tion frequency vary anywhere between 500 and 1500 Hz instead 
of one of the two previously used fixed frequencies (see 
Fig. 20) . 

2 . Differential Omega -Correction Information 

Goodman proposed use of an audio tone which varied 
in proportion to the Omega correction in the following 
manner ; 



Differential 
Omega Correction 


Audio 


Tone 


+50 


cec 


100 


Hz 


+ 25 


cec 


350 


Hz 


0 


cec 


600 


Hz 


-25 


cec 


850 


Hz 


-50 


cec 


1100 


Hz 



The audio tone variation proposed in this thesis is 
as follows : 



f. 
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Differential 
Omega Correction 


Audio Tone 


+50 


cec 


1500 Hz 


+25 


cec 


1250 Hz 


0 


cec 


1000 Hz 


-25 


cec 


750 Hz 


-50 


cec 


500 Hz 



The reasons for the changes are: 

a. Since there is an extremely stable 1000 Hz 
source available — an output at the back of the Tracor 599R 
Omega receiver — this can provide an extremely accurate source 
for a mid-scale calibration tone whose use will be described 
later. 

b. A minor change, but one that would be signifi- 
cant to a user, is to reverse the meaning of the high and 
low frequencies. It is natural for a user to think of a 
pointer deflection to the right of center as positive and to 
the left of center as negative. If Goodman's frequency 
schedule were fed into a frequency meter, mid-scale would 
correspond to zero correction, as in the newly proposed 
system, but high scale would correspond to a negative cor- 
rection; the reverse of the natural tendency of most users. 

3 . Coding Message Format 

Since the vehicle and modulation technique for 
transmitting the differential Omega information had been 
selected, it was necessary to decide in what format the 
information was to be transmitted. Both Goodman [7] and 
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Figure 21. Goodman's Differential Omega Message. 
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McKaughan [13] decided that the simplest, most easily decoded 
message was a correction for each of the four best hyperbolic 
LOPs through the monitor site (best as far as fix geometry 
and signal strength are concerned) . Goodman decided to 
transmit four differential Omega correction signals and five 
calibration signals in a one minute segment, repeating every 
five minutes (see Fig. 21) . The reason for chosing the one 
minute segment length was in order not to disrupt the normal 
radiobeacon chain transmission sequence. One radiobeacon in 
a net of five along the coast normally transmits during one 
minute out of five. The calibration tones were for users 
with an analog frequency meter which might require calibra- 
tion every 30 minutes. 

McKaughan 's coding message was slightly different 
(see Fig. 22). He decided to transmit the differential 
corrections continuously, repeating every minute for 29 
minutes, then transmitting a calibration tone for one minute. 
This would allow the alignment of analog frequency meters 
once every half-hour. McKaughan 's reason for changing from 
once every five minutes to continuous transmissions was 
given as the rapid change in Omega readings possible if a 
solar flare occurs at the end of the once per five minute 
transmission. For example [13] 

160° phase shift/40 min. = 4°/min. 

4° = 1.1 cec 

1 cec = 160 yds thus 1.1 cec = 176 yds error/min. 

5.5 cec = 880 yds error/5 min. 
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Figure 22. McKaughan's Differential Omega Message. 



Using a five minute old set of corrections could result in 
an error of 880 yds. This is clearly an unacceptable system 
error. 

Rapid changes of phase occur normally twice a day 
during the day-night, night-day transition period as can be 
seen in Figure 4. These changes can be as much as 

20 cec/hr = 1.66 cec error/5 min = 266 yds error/5 min 

(from Figure 4) 

Data taken for a typical LOP (Trinidad-Hawaii ) from skywave 
correction tables for the geographic position of the Naval 
Postgraduate School yielded the following results: 



24 hour average Acec/5 min = 0.73 cec error/5 min 
1-15 April 1972 

= 117 yds error/5 min 

largest Acec/5 min out of 24 hours = 2.3 cec error/5min 
10-1100 GMT (sunrise at Trin.) 

= 36 8 yds error/5min 

It can be seen that there is almost half as much incremental 
phase shift occurring during a five minute segment of the 
night to day transition period as can occur during a solar 
flare during the day. Therefore it was found necessary to 
transmit the differential Omega correction information on a 
once a minute schedule, as McKaughan had concluded, in order 
to derive the maximum benefit from its use during periods of 
normal, rapid phase change. In fact a navigator would prob- 
ably want correction information more often than once every 
five minutes if he were approaching a harbor entrance. 
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The differential Omega correction format proposed in 
this thesis can be seen in Figure 23. The only differences 
between this format and McKaughan's was the inclusion of a 
calibration tone every minute vice every thirty minutes and 
the method of identification of the correction lines. The 
inclusion of a calibration tone was deemed necessary since 
it is possible that a user might first decide to use the 
differential correction information immediately after the 
calibration segment was transmitted and would have to wait 
thirty minutes before he could be sure his frequency meter 
were properly calibrated. The identification as to which 
LOP the particular corrections applied would be published 
on the charts along with the symbol for the radiobeacon 
station and its operating frequency. For example; 

Radiobeacon 310 kHz 

'symbol LOP 1 B-D 

— LOP 2 A-D 

LOP 3 B-C 

LOP 4 A-C 

C. IMPLEMENTATION OF RADIOBEACON MODULATOR 
1 . Message Format Commutator 

The differential Omega message format described 
above requires a device that can shift any one of six dif- 
ferent audio tones to the transmitter input J102 (see Fig. 
24) . These tones are the 1000 Hz calibration and identifi- 
cation tone, LOPs one through four (audio tones of 500- 
1500 Hz) and no tone for the pause. A 100 Hz signal from 
the back of the Tracer 599R Omega receiver is put through • 
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Figure 23. Proposed Differential Omega Message Format 



two, divide by ten, counters so that a 1 Hz signal is then 

fed to the binary counters. The binary counter has the 

capability of counting up to 4095 seconds. Since the read 

only memory (ROM) used in the commutator can only accept 

the binary equivalent of 1023 different binary numbers at 

its input terminals , only the first ten outputs of the 

binary counter were used, A through . 

o y 

The ROM is a device that causes any combination of 

its four outputs, 01 through 04, to be activated by any 

combination of its inputs, A through A_ . A sample of the 

o y 

ROM programming is shown in Figure 25. ROM binary output 1, 
0001, corresponds to LOP correction 1, binary output 2,0010, 
to LOP correction 2 and so on. Binary output 5,0101, cor- 
responds to the 1000 Hz tone used for calibration and iden- 
tification and binary 0,0000, corresponds to ho tone for 
the pause. Since the ROM can accept up to 1023 binary 
inputs, this is equivalent to 1023 seconds or 17 min 3 
seconds. This means that the ROM can operate for seventeen 
one minute cycles before any of the input's numbers repeats 
itself. The three seconds left over are used to reset the 
binary counters. Count 1021 and 1022 will be the same as a 
pause and 1023 will be binary 6,0110. This will reset the 
binary counters to zero and allow the process to repeat 
again for the next 17 minutes. For the complete ROM truth 
table see Appendix A. 

The output of the ROM must be converted from binary 
to decimal numbers . The decimals numbers then correspond to 
the following list; 
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Figure 24. Block Diagram of Commutator Circuit. 
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Figure 25. Sample from Read Only Memory Truth Table. 
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ROM Output 
Binary Nuinber 



Function 



Decimal 

Equivalent 



0001 


1 


LOP 


#1 


0010 


2 


LOP 


#2 


0011 


3 


LOP 


#3 


0100 


4 


LOP 


#4 


0101 


5 


1000 Hz 
calibration 


0110 


6 


Reset 


0000 


0 


Pause 



From the binary to decimal converter all the pulses 
go to inverter circuits except number 6. The reason for 
this is that the outputs of the binary to decimal converter 
are all at plus 5 volts until that particular output is 
selected by the ROM, then it drops to zero volts. The 
inverter changes the outputs to zero volts until selected, 
then the output goes to 5 volts. The reset pulse doesn't 
go through an inverter since the reset pins of the binary 
counters require 5 volts there at all times. When the 5 
volts is removed then the counters reset. Since the six 
output of the binary to decimal converter is at 5 volts 
until selected, then this output can be connected to the 
reset pins of the binary counters and will cause them to 
reset when output 6 is selected by the ROM. 

After going through the inverter circuit the next 
step is the level translator to increase the voltage of the 
signal enough to operate the FET (Field Effect Transistor) 
switches. Each of the audio signals generated in the Omega 
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LOP corrector (described later) is fed into its own FET 
switch and all the outputs of the FETs are connected together 
and go to J102, the input jack of the audio modulation 
circuit of the transmitter. When the proper LOP has been 
selected by the ROM and decoded, then the switch correspond- 
ing to the LOP selected is allowed to feed through the FET 
to J102. All other FETs appear as open circuits to their 
signals. If no switch is selected then a pause occurs, 
i.e., there will be no audio modulation of the carrier. 

This portion of the circuit is the equivalent of five 
switches in parallel. This completes the commutator cir- 
cuitry. The complete schematic for the commutator is con- 
tained in Appendix B. 

2 . LOP Differential Correction Circuit 

As first envisioned, the LOP correction circuit was 
as shown in Figure 26. The Omega receiver output was a DC 
voltage that was directly proportional to the phase in cec 
of the LOP to which the receiver was tuned. This voltage 
was fed to a subtr actor circuit. A second signal provided 
by a stable DC source was also fed to the subtractor 
circuit. This was the reference to which the Omega signal 
was to be compared. The difference between the reference 
and the Omega receiver was an indication of the error of the 
Omega LOP from the correct value. This error voltage was 
then fed to a voltage controlled oscillator (VCO) which was 
used to generate an audio tone corresponding to the correc- 
tion signal. 
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a. Voltage Controlled Oscillator (VCO) 

Before any construction could begin on the over- 
all device, the frequency vs voltage characteristics of the 
VCO had to be obtained. The schematic diagram for the VCO 
is as shown in Figure 27. Using the formula 



f 

o 



2(V^ - V^) 



RlCiV 



(44) 



and a value of V^ of 12 volts, V =10.5 volts, R, = 4K ohms 

O -L 

Cj^ was determined to be 0.0625 yf. 

The circuit was built and tested and the results 
are shown in Figure 28. The linearity of frequency versus 
voltage was extraordinary. Frequency was measured to the 
nearest cycle and V^ , the control voltage, to the nearest 
millivolt. There were no observed non-linear regions in the 
graph within the designed range of 500-1500 Hz. 
b. Subtractor Circuit 

After determining the characteristics of the 
VCO, the next circuit attacked was the subtracter circuit. 
Using the circuit shown in Figure 29, the device was built 
and tested and found to have an offset voltage error. At 
the same time that the voltage offset error problem was 
being solved, it was realized that subtraction of the Omega 
signal from the reference voltage and its application to the 
VCO would cause the VCO to generate an audio frequency pro- 
portional to the error not the correction as desired. As 
an example assume that the voltage at the control terminal • 
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Figure 26. Block Diagram of LOP Correction Generator. 




Figure 27. Circuit Diagram for VCO. 
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Figure 28. Frequency Versus Voltage of VCO. 
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Figure 29. Subtracter Circuit. 




Figure 30. Block Diagram of Differential Omega 
Corrector Using Adder. 
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of the VCO, V / must be equal to 10.5 volts for the VCO to 
c 

generate 1000 Hz. Also assume that the voltage from the 
Omega receiver varies linearly from 0 to 1.0 volts corres- 
ponding to 0 cec and 100 cec respectively. To determine 
what voltage must be used for the DC reference when using a 
subtracter circuit, 

~ ^REF " ^OMEGA 

,^REF “ ^c ^OMEGA 

Assume that the correct LOP through the monitor site is 
50 cec; this corresponds to 0.5 volts. Thus using equation 
(46) 

= 10.5 + 0.5 = 11.0 volts 



So if the Omega receiver is actually receiving 50 cec then 
the VCO generates 1000 Hz as required. Now assume that the 
Omega reading rises to 60 cec or = 0.6 volts. Using 

equation (45) 

= 11.0 - 0.6 = 10.4 volts 
c 

As can be seen from Figure 28, a decrease in voltage means 
an increase in frequency. Assume that a of 10.4 volts 
generates 1100 Hz. This corresponds to a correction for the 
Omega LOP of +10 cec. This is an incorrect result since if 
the Omega receiver shows 60 cec it is 10 cec above the 
correct value and therefore the correction is -10 cec not 
+10 cec as generated using the subtracter circuit. 
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If an adding circuit replaced the subtracter circuit 



in Figure 26, then the applicable equations are 



V = V + V 
c REF OMEGA 



(47) 



V = V - V 
REF c OMEGA 



(48) 



Using the same assumptions as before (correct LOP 50 cec) , 
and equation (48) 



= 10.5 - 0.5 = 10.0 volts 

Kilir 



Using equation (47) 



V =10.0 +0.5 =10.5 volts 
c 



This is the correct voltage to generate 1000 Hz as 
required when the received Omega value is equal to the cor- 
rect value. Assume again that the Omega reading rises to 
60 cec or = 0.6 volts. Using equation (47)- 

V =10.0 +0.6 =10.6 volts 
c 

As before an increase in voltage corresponds to a decrease 
in VCO frequency. Assuming a of 10.6 volts equals 900 Hz 
then this corresponds to -10 cec, the proper magnitude and 
direction of a correction for an Omega reading which is 10 
cec too high. 

From these results the circuit in Figure 30 was 
devised. Upon careful consideration of this newly proposed 
circuit, it was seen that the same effect could be obtained 
if were a voltage floating with respect to the ground 

of the Omega receiver. By floating the negative terminal 
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of the reference voltage it could be placed at the positive 

terminal of the Omega receiver output voltage. This v/ould 

be the equivalent of placing two batteries in series, thus 

forming an adder circuit. By using a floating voltage it 

was possible to eliminate the two inverters and the inverting 

adder from the circuit in Figure 30 thereby eliminating 

three possible sources of nonlinearities. The chances for 

errors occurring in the summing circuit were high since the 

output would be close to 11 volts — definitely not a small 

scale signal and probably approaching the nonlinear portion 

of the operational amplifier's characteristics. 

c. Voltage Reference Circuit 

Since the subtracter and adder circuits had been 

eliminated, the next step was the design of the voltage 

regulator circuit. It was required that the voltage sources 

be extremely stable. The voltage of V , the VCO control 

c 

voltage, ranged from 10.325 to 11.811 volts for a frequency 
change of from 1500 to 500 Hz respectively. This range is 
only 1.486 volts or 1.486 x 10 ^ volts/cycle. Since accuracy 
to the nearest cycle of audio correction frequency was 
desired, this meant that the voltage references had to have 
a stability of better than 1.486 mv. The voltage regulator 
selected claims to have a long term stability of 0.1%/1000 hr 
or 10 mv/1000 hrs if = 10.0 volts. 1000 hours is over 

one month so this specification is not too important since 
the correction devices will at manned radiobeacon stations 
and will be checked three or four times daily by watchstanders 
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Figure 31. Schematic of Voltage Regulator Circuit. 
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Figure 32. Example of Cyclic Nature of Omega Corrections. 
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£or proper operation. Another parameter of the selected 
regulator is a temperature coefficient of 0.002%/°C or 
0.2 mv/®C assuming ~ 10.0 volts. The radiobeacons are 

generally housed in manned building and therefore the tem- 
perature would probably show only a few degrees change 
thereby eliminating this parameter as a source of error. 

In order to provide a floating voltage it was 
necessary to use a separate regulated power supply floating 
with respect to building ground to provide the 18.0 volt v"*" 
for the circuit shown in Fig. 31. The values of and R 2 
are calculated from the formulas 



V 



R-. = 



R.^ = 



out 



V ^ 

^ ref 

2000 V 



(2000) ohms 



out 



V . - V ^ 

out ref 



ohms 



(49) 

(50) 



is typically 1.58 to 1.68 volts and must be measured 
for each circuit to obtain the proper for the actual 

circuit on hand. 

If a particular Omega LOP passing through the 
monitor site were exactly 50 cec only one voltage regulator 
circuit would be needed. If the Omega LOP were any other 
value than 50 cec, assume 25 cec, then two voltage references 
are required because of the following: 



Omega LOP 
Received 


LOP Correction 
Required 


Frequency Generated 
Using One Voltage 
Regulator 


25 cec 


0 


1000 Hz 


50 cec 


-25 cec 


750 Hz 


75 cec 


-50 cec 


500 Hz 


15 cec 


+10 cec 


1100 Hz 


5 cec 


+20 cec 


1200 Hz 


85 cec 


+40 cec 


400 Hz (Incorrect) 



The last value in the list 400 Hz is incorrect. By examin- 
ation of Figure 32 it can be seen that 85 cec corresponds to 
an error of -40 cec or +60 cec. Since +50 cec was selected 
as the bounds for the LCP error by Goodman [7] , it was also 
used as one of the basic premises of this thesis. As can be 
seen by the above example and Figure 32 any reading greater 
than 75 cec would result in an incorrect correction frequency 
being generated for an assumed LOP of 25 cec if only one 
voltage reference were used. The general equations to deter- 
mine the upper and lower limits of the Omega readings from 
the correct values, assuming a +50 cec excursion are: 

Correct LOP less than 50 cec — 

Correct LOP + 50 cec = upper limit (51) 

0 cec = lower limit 

Correct LOP greater than 50 cec — 

Correct LOP - 50 cec = lower limit (52) 

100 cec = upper limit 
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The Omega receiver circuitry itself takes care of switching 
from 100 to 0 cec when phase increases through 100 cec and 
0 to 100 cec when phase decreases through 0 cec. 

There must be some means of switching another 
voltage source into the circuit when exceeds the 

limit set by equation (51) or (52) . The switching can be 
accomplished by a voltage comparator using the circuit shown 
in Figure 33. A comparator circuit is basically an open 
loop operational amplifier with no degenerative feedback and 
a gain of perhaps as much as 10^. A difference between 

^OMEGA ^COMPREF millivolt would be enough to 

drive the output of the comparator into saturation (+V^^) . 



This causes very sharp switching when V 



OMEGA 



reaches 



^COMPREF gradual output voltage change at the com- 



parator as approaches When 



IS 



greater than the output of the comparator is +V . 

UJyLCiUA cc 

''COMPREF *^*'®*' ''omega "''cc' 

The output of the comparator is fed to a voltage 

follower in order to isolate the comparator circuit from the 

heavy drain caused by the relay. The input resistance of 

the voltage follower is nominally 400 MEGohms while its 

output resistance is much less than one ohm. The output 

voltage of the voltage follower is an exact duplicate of its 

input voltage hence the name voltage follower. When the 

comparator — voltage follower output is +V , meaning that 

^COMPREF greater than then the diode is forward 

biased and the relay is closed allowing to be used in 
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Figure 33. 



the VCO circuit. If is greater than then 

the output of the comparator — voltage follower is so 

the diode is reverse biased and the relay opens connecting 
Vref 2 to the VCO circuit. The capacitor in parallel with 
the relay is to eliminate relay chatter. 

As an example to help clarify the need for two 
voltages reference assume that: 

(a) of the VCO must be 10.5 volts to generate 
1000 Hz 

(b) V of 10 volts = 1500 Hz, V of 11.0 volts = 

c c 

500 Hz 

(c) 0 cec = of 0 volts, 100 oec = 

of 1.0 volts 

(d) The correct value of the LOP of interest is 

25 cec = 0.25vDlts.) 

Then using equation (48) to solve for V„„„ 

Kljr 

^REF ^ ^c " ^OMEGA (48) 

Vrefi = 10.5 - 0.25 volts 



Vrefi = 10.25 volts 



Equation (51) must now be used to determine the value of 
the upper limit of V,_..„„_ above which switching must occur 

UlXlri LjA 

Correct LOP + 50 cec = upper limit (51) 

25 cec + 50 cec = 75 cec 



75 cec corresponds to 0.75 volts, thus = 

COMP REF 

0.75 volts. A value of 75 cec could require a correction 
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of either +50 cec or -50 cec since it is the end point 
either way from 25 cec as shown in Figure 32. Therefore to 



determine the value of given that 75 cec equals a 

+50 cec correction (forcing V to equal 10.0 volts so that 

c 

the VCO will generate 1500 Hz or +50 cec) ; 



^REF2 

^REF2 

^REF2 



V - V 
c OMEGA 

10.0 - 0.75 volts 
9.25 volts 



Now if is less than 0.75 volts then V_„„, 

OMEGA REFl 

equals 10.25 volts will be used as the voltage reference for 
the VCO and if is greater than 0.75 volts, 

equals 9.25 volts will be used as the reference. 

d. Final Differential Omega Correction Circuit 
Figure 34 shows the LOP correction circuit in 
its final form. The Omega signal, is first routed 

to the voltage comparator circuit where it is 
compared to ^^OMPREF which of two voltage sources, 

^REFl '^REF2' ^iii switched by the relay to the adder 
circuit. The adder circuit merely consists of the floating 
reference voltages, or Vj ^^2 being placed in series 

with the grounded and their combined sum then being 

fed as to the VCO. The output of the VCO goes through a 
voltage follower circuit to buffer the output from the 
effects of the transmitter. This was. found necessary because 
when the VCO was connected to the receiver, the frequency of 
the VCO changed 20 Hz. The VCO has its own voltage regulator 
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circuit providing it is a of 12 volts. The VCO requires 
its own power supply because, as can be seen from solving 
equation (44) , the output frequency is extremely dependent 
on the value of v'*'. Solving equation (44) for the differ- 
ential change in f^ versus v"*" gives: 





1000 

12 volts 



0.58 Hz 
millivolt 



Only a 2 millivolt change in results in a 1.06 Hz change 
in output frequency or 0.1 cec when converted to LOP correc- 
tion terms. For the complete schematic diagram of the dif- 
ferential Omega LOP corrector see Appendix C. 
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V. TESTING AND ANALYSIS 



A. RESULTS OF TESTING WITHOUT TRANSMITTER 

When construction of the differential Omega LOP corrector 
was completed the next step was a thorough testing of the 

r 

device. It was decided to test initially the device without 
using the radiobeacon transmitter. This would give an 
indication of the accuracy of the corrector without errors 
from interfacing it with the transmitter. The first tests 
were conducted using the equipment as arranged in Figure 35. 
The output of the Omega receiver was recorded on the red 
channel of a dual channel strip chart recorder moving at 
0.1 inches per minute. The chart paper was ten inches wide 
and divided into 1 inch major and 0.1 inch minor divisions. 
Zero scale corresponded to zero cec and full scale corres- 
ponded to 100 ceC/ thus each minor division was 1 cec. The 
Omega receiver output was also fed to the differential Omega 
LOP corrector where it was used to generate the audio cor- 
rection frequency. The audio correction frequency was fed 
to a Hewlett Packard 5210A analog frequency meter. This 
meter has a recorder output which was fed to the blue channel 
of the chart recorder. Zero scale on the blue channel cor- 
responded to 500 Hz and full scale corresponded to 1500 Hz. 
Thus each minor division corresponded to 10 Hz. 

Figure 36 shows a sample of a chart recording for 2 March 
1972. It graphically illustrates the magnitude of phase 
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Figure 35. First Differential Omega LOP Corrector Test 
Layout Without Transmitter. 
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Figure 36. Effect of Solar Flare on Phase of Omega Signal and Differential 
Omega LOP Corrector, 2 March 1972. 



change due to a small anomalous disturbance. There is a 
sudden phase decrease of approximately 10 cec in five minutes, 
probably due to a solar flare. The lower curve is the out- 
put of the differential Omega LOP corrector as recorded from 
the analog frequency meter. This curve shows an increase in 
frequency of 100 Hz which is the change in frequency required 
to correct a 10 cec decrease in phase. 

The data was normally analyzed once a day. The charts 
were read every ten minutes and the received Omega LOP was 
recorded as well as the frequency being generated by the 
corrector at that time. Both recorded lines were interpolated 
to the nearest 0.1 of a minor division or 0.01 inch. An 
Omega LOP or frequency measurement read to the nearest 0.01 
inch resulted in accuracies to the nearest 0.1 cec for the 
former and 1 Hz for the latter. These two numbers were then 
placed on computer data cards and run through the Naval Post- 
graduate School's IBM 360 computer. The program used is 
located in Appendix D. 

The computer program calculated the difference between 
the Omega LOP being received and the correct value. Based 
on this difference, it then computed what required frequency 
the differential Omega LOP corrector should have been gener- 
ating to correct the received LOP to the proper value. After 
this the computer read the frequency actually generated from 
the data card and computed the difference between what the 
LOP corrector generated and what it should have generated. 

This difference resulted in a frequency error which was 
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converted to cec error by dividing by ten (A frequency 
1000 Hz = Acec 100) . The average value of the cec error was 
obtained using equation (53) 



Average error cec 



n 

I, (errcec) . 

i^l i 

n 



(53) 



The standard deviation of the cec error was obtained using 
equation (54) . 



o 

error cec 



/ n 

nl Z (errcec) 
\i=l 




/ "" 

I r (errcec) . 

\i=l ]■ 




(54) 



2 



n 



where errcec^ is the individual error of each data point in 
cec and n is the number of observations made. A sample com- 
puter printout is contained in Appendix E. 

The assumed Omega LOP through the antenna location used 
for this experiment was obtained from Mr. E. R. Swanson of 
NELC, San Diego. The Trinidad-Hawaii (B-C) 10.2 kHz LOP for 
this location is 8.5 cec. 

Figure 37 is a summary of the results of the computer 
output. It contains results for and need 



for two reference voltages was explained earlier. V 



REFl 



was 



used for received Omega signals of less than 58.5 cec. Above 

58.5 cec was switched into the circuit instead of 

REF2 

^REFl* 

In the case of non-correlated lines of position, the 
standard deviation of the fix error is usually taken as the 
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Figure 37. Summary of Accuracy of. Differential Omega Corrector Using 
Strip Chart Recorder without Transmitter. 



square root of the sum of the squares of the two LOP errors 
divided by the sine of the crossing angle between the two 
LOPs [19]. The crossing angle was shown earlier to never be 
less than 60®, and 72° may be assumed as a typical value. 
Assuming both LOPs can be corrected to the same accuracy, 
then the standard deviation of the fix can be determined 
from equation (56) . 

a /T" 

^ errcec 

°fix error sin 72° ^ ' 



fix error 



= (o 



errcec 



) (1.49) 



(56) 



Equation (56) gives the fix error in cec which can be con- 
verted to yards by multiplying by 160 (Icec = 480 ft. =160 yds 
on the baseline) . The ^f^xerror yards for and 

51 yards for 

After 16 March 1972, a new data acquisition technique 
was employed (see Fig. 38). The analog frequency meter 
was replaced by a Hewlett Packard digital frequency meter. 

The Omega receiver voltage, measured by a Hewlett 

Packard digital voltmeter. Instead of recording the data 
on an analog strip chart, the data was then recorded on 
three inch wide adding machine paper by a Hewlett Packard 
560A digital recorder. The paper printout was not in binary 
numbers but ordinary decimal numbers. The frequency was 
recorded to the nearest Hz and to the nearest 0.10 cec. 

The paper tape information was then transcribed to data cards 
and the same computer program above was used to analyze the 
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data. Data points were obtained every six minutes by a 
clock operated cam. This cam fed a signal to the frequency 
meter causing it to sample the frequency. When the fre- 
quency meter completed its one second count, it commanded 
the printer to print the frequency and the digital volt- 
meter reading. Figure 38 contains the results of the com- 
puter runs. During the period covered by Figure 39, 16-25 

March 1972, 0 _. of was 21 yards and was 

fixerror REFl ^ REF2 

30 yards. The reduction of *^fixerror Figure 37 to 

Figure 39 was attributed to the fact that the data inter- 
polation was being done automatically instead of manually. 

B. RESULTS OF TESTING WITH TRANSMITTER 

After sufficient data had been collected to determine 
what the accuracy of the differential Omega LOP corrector 
was by itself, it then had to be tested as a transmitter 
modulator (see Fig. 40). The output of the corrector was 
fed to J102, the audio m.odulation input jack of the 
T-747A/FRN radiobeacon transmitter. The transmitter was 
connected to a dummy antenna. A wire antenna for the 
AN/WRR-3A low frequency receiver was placed near the dummy 
antenna and the stray emissions were picked off and used 
as receiver input signals. The signals were demodulated and 
the audio output of the receiver was fed to the digital 
frequency meter. The same cam timer as before resulted in 
data points being recorded every six minutes. The output , 
of the printer was typed onto data cards and run through 
the computer. The results are shown in Figure 41. 
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Figure 38. Second Differential Omega LOP Corrector Test 
Layout Without Transmitter. 
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Figure 39. Suinmary of Differential Omega Corrector Accuracy 
Using Digital Recorder Without Transmitter. 
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Figure 40. Differential Omega LOP Corrector Test Circuit with Transmitter. 
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41. Summary of Accuracy of Differential Omega Corrector with Transmitter. 



Unexpectedly, there was very little loss of information 



by using the transmitter-receiver combination. The 

o_c- was 24 yards and was 35 yards, 

fixerror REFl ^ REF2 ^ 

This was only a 4 yard (0.2 cec) increase over operation 



without the transmitter for V 



REFl 



and a 5 yard (0.03 cec) 



increase for V 



REF2* 



C. ANALYSIS OF RESULTS 



1 . Test Results 

The circular probable error (CPE) is considered to 
be the radius of a circle, centered at the mean position of 
a group of measurements, that encloses 50 percent of, the 
measurements . 



CPE 



1.18 a 

errcec 

sine crossing angle 



(57) 



Ninety-five percent of the fixes will lie within a circle 
of radius: 



2.45 0 



95% circular error = 



errcec 



sine crossing angle 



(58) 



Using equation (57) and (58) and the averages of 
Figure 41 the results are: 



^REFl = = 0.12 cec = 19 yards 

CPE = '"^ '' ^0 ~ 95* “ 0.19 cec = 30 yards 



Average CPE = 25 yards 
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95% Circular error V. 



2.45 (0.1) 
REFl 0.95 



0.25 cec = 41 yards 



95% Circular error V. 



REF2 



2.45 (0.15) 
0.95 



=0.39 cec = 62 yards 



Average 95% Circular error = 52 yards 

These results mean that at the monitor site, an 
erroneous Omega fix can be corrected to within a circle of 
52 yards radius 95% of the time and a circle of 25 yards 
radius 50% of the time. This is a tremendous improvement 
over the generally conceded 1 to 2 mile error of ordinary 
Omega using skywave correction tables. 

2 . Digitizing Errors 



some errors arise. Accuracies of digital devices are 
generally held to be + one in the last digit read. The 
standard deviation of a digital device assuming that the 
likelihoods that the device will read one high, one low or 
right on are equal gives : 



Anytime a digital measurement is made inevitably 




(59) 



where n = 3 in this case 




i 



a = (0.75)^ 



a = 0.866 counts 



a 



cec 



0.0866 or 0.09 cec 
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This means that 0.09 cec or 14 yards of the error 
for one LOP attributable to the differential Omega corrector 
could be due to digitizing error and not error in the cor- 
rector device itself. This would mean that the digitizing 
error by itself could cause a CPE of 0.11 cec or 18 yds and 
a 95% circular error 0.23 cec or 37 yds. It is felt that 
while there is some digitizing error in the results pre- 
sented in Figure 41, 0.09 cec is too large a value. The 
probability of being one count high is not necessarily the 
same as being one count low, and neither is the same as 
the probability of being right on but in general each is 
smaller. This would reduce the standard deviation of the 
digitizing error to something smaller than 0.09 cec but it 
is an undeterminable amount since the exact probabilities 
of each of the digitizing errors cannot be obtained. 

3. Receiver Errors 

All the results discussed previously have been using 
the assumption Ihat the receiver at the monitor site and the 
mobile receiver would read the same if they were placed side 
by side and fed from the same antenna. This however is not 
the case for the receivers used in this experiment, the 
Tracer 599R. NELC conducted tests on three 599R receivers 
and concluded that they could be expected to agree to within 
1 cec about half the time; to within 2 cec most of the time 
and to within 3 cec essentially all the time. These are not 
statistically analyzed results but are results of an intui- 
tive analysis by the author of Ref. 22 [22] . Based on this 
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information, if Tracor 599R receivers were used as the 
monitor receiver and the mobile receiver, and both were 
receiving the same station and were theoretically on the 
same phase difference, the LOPs obtained would agree within 
160 yds half the time, 320 yards most of the time and 
within 480 yds essentially all the time. This leads one to 
believe that in any experiment using two receivers, extensive 
analysis of their measurement errors should be undertaken 
first. If two receivers had been used in this experiment 
the 95% circular error of 52 yards would have been lost in 
the errors caused by the disparities between the two 
receivers. Tests conducted on two Tracor 599R receivers at 
the Naval Postgraduate School showed that when both receivers 
were tracking the same signal an error of as much as one cec 
occurred (see Fig. 42). Additional tests to show the 
difference between two channels on the same receiver track- 
ing the same signal resulted in an error between the two 
channels of as much as 1.5 cec (see Fig. 43 ). 

If Tracor 599R Omega receivers were again used as 
monitor and mobile receivers and both were receiving the 
same station but were theoretically on different phase LOPs, 
there could be as much as 1.4 cec difference between the two 
readings. This error is caused by what is known as s-curve 
error. s-curve error is caused by the receiver not display- 
ing the correct phase for the phase fed into it. For 
example if a 0 cec signal were fed info the receiver it 
might read 0.3 cec. If a 50 cec signal were fed in, the 
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receiver might read 50.8 cec and so on. The reason for this 
error being called S-curve error is because the shape of the 
curve of error vs input is generally "S" shaped. S-curve 
error is different from the error described in the preceding 
paragraph because there the signals being compared were the 
same phase whereas in this paragraph the signals have dif- 
ferent phases. S-curve errors on the three NELC Tracer 
receivers ran from 0.6 to 1.4 cec peak to peak [22]. S-curve 
errors are believed to be caused by low level leakage of the 
internally generated 1 kHz signal in or around the tracking 
filters . 

It should be noted that the receiver errors described 
in this section cannot be removed by differential Omega. 

4 . LOP Errors Caused by Ship Velocity 

A ship moving at 15 knots along the baseline between 
two stations, X and Y is covering 8.34 yards/sec. The time 
difference between the stop of X's transmission and start 
of Y's could be over four seconds due to the Omega signal 
format shown in Figure 2. Therefore the ship could travel 
25.36 yards toward station Y in four seconds. Assume that 
if the ship were standing still and the phase of X = 50 cec 
and Y = 25 cec, the phase difference then, X-Y = 25 cec. 

Now assume the ship is moving and station X's phase was 
received as 50 cec and stored until Y comes on the air. 

When Y's transmission arrives, the phase will not be 25 cec 
as before but will be 24.84 cec because the vessel has moved 
25.36 (0.16 cec) toward Y during the four second interval 
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between X's and Y's transmissions. Thus X-Y = 25.16 cec 
causing a 25.36 yd. error in the LOP. This type of error 
can easily be removed by providing a biasing voltage for 
the readout circuit that can be adjusted to remove the 
error from each LOP. The bias voltage would be proportional 
to the velocity of the ship perpendicular to the LOPs in 
use and could be either positive or negative depending on 
whether the vessel was traveling toward or away from the 
Omega station. 

D. ADVANTAGES AND DISADVANTAGES OF USING THE COAST GUARD 

RADIOBEACON AS A DIFFERENTIAL OMEGA CORRECTION TRANSMITTER 

1. Advantages 

The primary reason for the implementation of the CG 

t 

radiobeacon differential Omega system is cost. The radio- 
beacon net is a fact of life', it is here and operating. The 
Omega system will be operational world-wide by the end of 
1974. The use of the radiobeacon as a means of transmitting 
corrections would allow one system. Omega, to be used for 
both mid-ocean and coastal navigation. 

As far as shipboard users are concerned, the only 
additional items needed for the minimum differential Omega 
installation would be a radiobeacon receiver and a frequency 
meter. Any ship using Omega as a navigation aid would 
undoubtably have a LF receiver aboard. The proposed system 
could be used by sophisticated computer oriented users, or 
the unsophisticated navigator without changing the message 
format . 
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The Raydist and Lorac systems both suffer from the 
fact that not only are expensive transmitter sites (three 
or four per chain) required every 200 n.mi. up and down the 
coast, but a continued operating cost would be incurred. 

The user would also be required to purchase special receivers 
in order to use these systems. In addition the range of 
these systems is reduced at night. 

The Sercel system and Micro-Omega also require 
additional transmitter sites (one per station) to be built — 
with their continued operating cost--plus special receiving 
units on the ships. No cost figures are available for the 
Micro-Omega system hardware. The Sercel differential 
correcting receiver and standard Omega receiver will cost 
approximately one-third more than the Sercel Omega receiver 
alone. The receiver will cost approximately $5000 without 
the differential correction receiver. 

Loran A staitons cost $250,000 each per year to 
operate. When the Omega system becomes operational with 
all stations at full power, it appears superfluous to con- 
tinue operating Loran A stations. Estimated operating 
costs for Loran A over the next 30 years is 360 million 
dollars for the U.S. alone [24]. Loran C operation costs 
are even higher and furthermore the receivers cost at 
least $20,000 for an automatic tracking receiver compared 
to $4000 to $7000 for an Omega receiver. 

2. Disadvantages 

One definite disadvantage common to all differential 
Omega systems is the basic lane ambiguity problem suffered 
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by the Omega system itself. The proposed radiobeacon system, 
the Sercel system and Micro-Omega are alike in that they 
cannot help solve the ambiguity problem. Raydist and Lorac 
are also subject to the lane ambiguity problems of Omega 
only on a smaller scale. Of course Loran A and C suffer no 
ambiguities so they have the advantage over Omega in this 
respect. 

The Sercel and Micro-Omega systems have the advantage 
of being inherently automatic and therefore simpler for the 
navigator to operate. Demands on the attention of the 
navigator are likely to be greatest when he needs the addi- 
tional accuracy of differential Omega and the simpler the 
system is to operate the better. Corrected position data 
should be available for convenient visual readout. For 
this reason an automatic system would be favored. This 
advantage for the Sercel and Micro-Omega systems would be 
eliminated in the automatic shipboard installation using 
the radiobeacon system described earlier. Whereas the 
Sercel and Micro-Omega systems are real time correction 
schemes, the radiobeacon LOP correction information could 
be as much as one minute old. 

The use of the proposed continuous differential 
Omega transmissions from a radiobeacon will disrupt the 
radiobeacon chain transmission technique. This will require 
the reallocation of the frequencies in the radiobeacon band 
to allow the differential Omega transmission station to 
transmit continuously. The frequency reallocation would 
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permit the remaining radiobeacon transmitters to operate as 
they are presently but on different frequencies than the 
differential Omega station. An additional cost necessary to 
implement the proposed system would be the purchase of a 
power amplifier for the differential Omega radiobeacon 
station. The purpose of this amplifier would be to boost 
the range of the radiobeacon from its present 125 n.mi. to 
200 n.mi. or more. 

As the radiobeacons transmit in the low frequency 
range, receivers are subjected to a great deal of atmos- 
pheric noise particularly during thunder storms. The ground 
wave is rapidly attenuated and night and skywave inter- 
ference can occur between the ground and skywaves . 

The disadvantages presented above however are far 
outweighed by the economic and operational advantages of 
implementing the system using the presently installed 
radiobeacons . 
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VI. SUMMARY 



The Omega navigation system will be implemented world- 
wide by the end of 1974. This system, using only eight 
transmitters, will be capable of providing navigators their 
choice of either 10 or 15 lines of position anywhere in the 
world. The Omega system will be able to secure world-wide 
navigation with only eight stations because the very low 
frequencies (10-14 kHz) employed show very little attenua- 
tion, good stability and excellent predictability. However 
VLF frequencies are subject to unpredictable (anomalous) 
effects caused primarily by the sun. Differential Omega is 
the process of observing the anomalous effects and dissem- 
inating correction information to users in the differential 
area. Persons using the differential information can 
theoretically expect to have fix accuracies of 0.26 to 0.5 
n.mi. at 200 n.mi. from the monitor site. The error 
decreases as the distance to the monitor site decreases. 
These errors are a great improvement over the generally 
expected 1 to 2 n.mi. fix errors associated with the use of 
ordinary Omega and published skywave correction tables. 

The presently operating system of U.S. Coast Guard 
radiobeacons was proposed by Goodman and McKaughan as a 
logical and inexpensive method of transmitting differential 
Omega information to users nearby. Their proposed systems 
were examined and improvements suggested. In addition, two 
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other differential Omega systems were examined and their 
advantages and disadvantages discussed. It was concluded 
that the radiobeacon system provided the most economic 
method of implementing a differential Omega system. 

Other coastal navigation systems were also examined but 
the use of differential Omega along with Omega appears to 
be the best solution since it can provide both mid-ocean 
and coastal navigation information with a minimum expendi- 
ture for. both the user and the government. The radiobeacon 
system has another great advantage in that the differential 
information can be extracted from the radiobeacon using as 
sophisticated a system as the user can afford — starting with 
only an analog frequency meter, through a more expensive 
digital frequency meter, then to a special purpose computer ' 
and on to a sub-routine on a large shipboard computer for 
the highly sophisticated user. 

A differential Omega LOP device was built and tested in 
a closed loop transmission system to determine its accuracy. 

A CPE of 25 yards and a 95% circular error of 52 yards were 
obtained at the monitor site. These results were analyzed 
and it was found that up to 18 yds of the CPE and 37 yds 
of the 95% circular error could be attributed to digitizing 
errors in the measuring equipment. 

Receiver errors are a major source of error in trying to 
optimize a differential Omega system. It was shown that 
receiver errors can be oh the same order of magnitude as the 
spatial correlation errors out to distances of about 100 n.mi. 
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and are a considerable percentage in size of the spatial 
correlation error out to much greater distances. The 
receiver errors are added to the errors predicted by the 
spatial correlation model causing gross errors compared to 
the errors of the differential Omega LOP corrector construc- 
ted for this thesis. For this reason it has been concluded 
that to derive the maximum benefit from the use of differ- 
ential Omega that a receiver with an S-curve error and a 
repeatability error of no larger than 0.1 cec will be 
required for use both at the monitor and mobile sites. 

Ship velocity errors can be removed by the addition to the 
receiver of LOP velocity correctors. 

It is strongly felt that any decision on a differential 
Omega system should be delayed until the entire 8 station 
Omega ^stem is operating at the rated 10 kW of power. Only 
then can an intelligent decision as to the full capabilities 
of Omega and differential Omega be made. To dismiss the 
many advantages of the use of differential Omega based on 
the presently available low power station would be a serious 
error. 

Future research efforts should include actual "on the 
air" testing in some location where there is an accurate 
navigation system available against which to measure the 
proposed system's error. This study should include both 
day and night testing to identify the errors caused by the 
interference caused by skywaves as well as low signal to 
noise ratios for the radiobeacon signal. 
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The system proposed in this thesis limits the correction 
Information to only +50 cec. Fifty centicycles is not 
adequate correction in all cases. For example, from 
0030-0415 GMT, 1-15 April 1972 the skywave correction is 
greater than -50cec for the geographic position of the 
Naval Postgraduate School using Trinidad-Hawaii (B-C) at 
10.2 kHz. If the error in the LOP exceeds 50 cec then the 
sign of the correction message abruptly changes as the 
increasing error goes through 50 cec. This can cause an 
error of one complete lane for a vessel which moments ago 
had its ambiguity problem solved. Future researchers should 
consider some means of ' identifying these occurrences. 

The correction message is presently restricted to one 
minute lengths. This was done originally in order to allow 
the differential Omega transmission to fit into the one 
minute transmission setments presently being used by the 
five transmitter radiobeacon chains. However this restric- 
tion has since been removed because it was realized that 
differential Omega transmissions spaced five minutes apart 
were not sufficiently accurate. The correction message 
should now be examined by future researchers from the stand- 
point of the best format from the users point of view. Any 
changes made to the format proposed in this thesis should 
strike a balance between user acceptability and minimum time 
between correction messages for any particular LOP. The 
minimum time criterion is in order to reduce the errors 
arising during normal, rapid phase changes during the 
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day-night, night-day transition period. Any new message 
format can easily be incorporated into the commutator circuit 
designed in thesis by changing the read only memory truth 
table. 

One last area of future research could be the design 
of the special purpose computer to be used in the automatic 
differential Omega correction system on board larger ships . 
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